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Abstract 

The current standard theory of the origin of the Moon is that the 
Earth was hit by a giant impactor the size of Mars causing ejection 
of iron poor impactor mantle debris that coalesced to form the Moon. 
But where did this Mars-sized impactor come from? Isotopic evi- 
dence suggests that it came from lAU radius in the solar nebula and 
computer simulations are consistent with it approaching Earth on a 
zero-energy parabolic trajectory. But how could such a large object 
form in the disk of planetesimals at lAU without colliding with the 
Earth early-on before having a chance to grow large or before its or 
the Earth's iron core had formed? We propose that the giant impactor 
could have formed in a stable orbit among debris at the Earth's La- 
grange point L4 (or L5). We show such a configuration is stable, even 
for a Mars-sized impactor. It could grow gradually by accretion at L4 
(or L5), but eventually gravitational interactions with other growing 
planetesimals could kick it out into a chaotic creeping orbit which we 
show would likely cause it to hit the Earth on a zero-energy parabolic 
trajectory. This paper argues that this scenario is possible and should 
be further studied. 



1 Introduction 

The currently favored theory for the formation of the Moon is the giant 
impactor theory formulated by Hartmann & Davis (1975) and Cameron & 
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Ward (1976). Computer simulations show that a Mars-sized giant impactor 
could have hit the Earth on a zero-energy parabolic trajectory, ejecting im- 
pactor mantle debris that coalesced to form the Moon. Further studies of 
this theory include (Benz, Slattery, & Cameron 1986, 1987; Benz, Cameron, 
& Melosh 1989; Cameron & Benz 1991; Canup & Asphaug 2001; Cameron 
2001, Canup 2004; Stevenson 1987). We summarize evidence favoring this 
theory: (1) It explains the lack of a large iron core in the Moon. By the late 
time that the impact had taken place, the iron in the Earth and the giant 
impactor had already sunk into their cores. So, when the Mars-sized giant 
impactor hit the Earth in a glancing blow, it expelled debris, poor in iron, 
primarily from mantle of the giant impactor which eventually coalesced to 
form the Moon (cf. Canup 2004, Canup2004B). Computer simulations (as- 
suming a zero-energy parabolic trajectory for the impactor) show that iron 
in the core of the giant impactor melts and ends up deposited in the Earth's 
core. (2) It explains the low (3.3 grams/cm^) density of the Moon relative 
to the Earth (5.5 grams/cm^), again due to the lack of iron in the Moon. 
(3) It explains why the Earth and the Moon have the same oxygen isotope 
abundance - the Earth and the giant impactor came from the same radius in 
the solar nebula. Meteorites originating from the parent bodies of Mars and 
Vesta, from different neighborhoods in the solar nebula have different oxygen 
isotope abundances. The impactor theory is able to explain the otherwise 
paradoxical similarity between the oxygen isotope abundance in the Earth 
combined with the difference in iron. This is perhaps its most persuasive 
point. (4) It explains, because it is due to a somewhat unusual event, why 
most planets (like Venus and Mars, Jupiter and Saturn) are singletons, with- 
out a large moon like the Earth. Competing ideas have not had comparable 
success. For example, the idea that the Earth and the Moon formed together 
as sister planets in the same neighborhood fails because it doesn't explain the 
difference in iron. Whereas the idea that the Moon formed elsewhere in the 
solar nebula and was captured into an orbit around the Earth fails because 
its oxygen isotope abundances would have to be different. That a rapidly 
spinning Earth could have spun off the Moon(from mantle material) is not 
supported by energy and angular momentum considerations, it is argued. 

Still, the giant attractor theory has some puzzling aspects. Planets are 
supposed to grow from planetesimals by accretion. How did an object so large 
as Mars, form in the solar nebula at exactly the same radial distance from 
the Sun without having collided with the Earth earlier, before it could have 
grown so large. Indeed, such must have been the case during the formation 
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of Venus and Mars for example. It's also hard to imagine an object as large 
as Mars forming in an eccentric Earth-crossing orbit. One might expect 
large objects forming in the solar nebula to naturally have nearly circular 
orbits in the ecliptic plane, like the Earth and Venus. Besides, a Mars- 
sized object in an eccentric orbit would not be expected to have identical 
oxygen abundances relative to the Earth, and would collide with the Earth 
on a hyperbolic trajectory not the parabolic trajectory that the successful 
computer simulations of the great impact theory have been using. (Recent 
collision simulations by Canup(2004) place an upper limit of 4 km/s for the 
impactor's velocity-at-infinity approaching the Earth, setting an upper limit 
on its eccentricity of ~ 0.13.) The Mars-sized object needs to form in a 
circular orbit of radius 1 AU in the solar nebula but curiously must have 
avoided collision with the Earth for long enough for its iron to have settled 
into its core. Is there such a place to form this Mars-sized object? 

Yes: the Earth's Lagrange point L4(orL5) which is at a radius of 1 AU 
from the Sun, with a circular orbit 60° behind the Earth(or 60° ahead of 
the Earth for L5). After the epoch of gaseous dissipation in the inner solar 
nebula has passed we are left with a thin disk of planetesimals interacting 
under gravity. The three-body problem shows us that the Lagrange point 
L4(or equivalently L5) for the Earth is stable for a body of negligible mass 
even though it is maximum in the effective potential. Thus, planetesimals 
can be trapped near L4 and as they are perturbed they will move in orbits 
that can remain near this location. This remains true as the Earth grows 
by accretion of small planetesimals. Therefore, over time, it might not be 
surprising to see a giant impactor growing up at L4(or L5). In the Discussion 
section we argue that there are difficulties in having the giant impactor come 
from a location different from L4(or L5). 

Examples of planetsimals rciiiciiiiins" cit Let grange points of other bodies 
include the well known Trojan asteroids at Jupiter's L4 and L5 points. As 
another example, asteroid 5261 Eureka has been discovered at Mars' L5 point. 
(There are five additional asteroids also thought to be Mars Trojans: 1998 
VF^i, 1999 C/J7, 2001 DH^T, 2001 FG2A, 2001 FRi27.) The Saturn system 
has several examples of bodies existing at the equilateral Lagrange points 
of several moons, which we discuss further in the "Note added" after the 
Discussion section. 

We propose that the Mars-sized giant impactor can form as part of de- 
bris at Earth's L4 Lagrange point. (It could equally well form at L5, but 
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as the situation is symmetric, we will simply refer in the rest of the paper, 
unless otherwise indicated, to the object forming at L4; the argument be- 
ing the same in both cases). As the object forms and gains mass at L4, we 
can demonstrate that its orbit about the Sun remains stable. Thus, it has 
a stable orbit about the Sun, and remaining at L4 keeps it from collision 
with the Earth as it grows. Furthermore, this orbit is at exactly the same 
radius in the solar nebula as the Earth so that its oxygen isotope abundances 
should be identical. It is allowed to gradually grow and there is time for its 
iron to settle into its core, and the same also happens with the Earth. The 
configuration is stable providing the mass of the Earth and the mass of the 
giant impactor are both below .0385 of the mass of the Sun, which is the 
case. But eventually, we numerically demonstrate that gravitational pertur- 
bations from other growing planetesimals can kick the giant impactor into 
a horseshoe orbit and finally into an orbit which is chaotically unstable in 
nature allowing escape from L4. The giant impactor can then enter an orbit 
about the Sun which is at an approximate radial distance of 1 AU, which 
will gradually creep toward the Earth; leading, with large probability, to a 
nearly zero-energy parabolic collision with the Earth. Once it has entered 
the chaotically unstable region about L4, a collision with the Earth is likely. 
We will discuss this phenomenon in detail in Section 4. For references on the 
formation of planetesimals and related issues, see (Goldreich 1973; Goldreich 
& Tremaine 1980; Ida & Makino 1993; Rafikov 2003; Wetherill 1989). We 
are considering instability of motion near L4 due to encounters by planetes- 
imals. (The instability of the Jupiter's outer Trojan asteroids due to the 
gravitational effects of Jupiter over time studied in Levinson, Shoemaker E. 
M. & Shoemaker C. S. (1997), is a different process.) 

Horseshoe orbits connected with the Earth exist. In fact, an asteroid 
with a 0.1 km diameter, 2002 AA2g, has recently been discovered in just 
such a horseshoe-type orbit which currently approaches the Earth to within 
a distance of only 3.6 million km (Conners et al. 2002). Horseshoe orbits 
about the Sun of this type are also called Earth co-orbiting trajectories, which 
are in 1:1 mean motion resonance. An interesting pair of objects in horseshoe 
orbits about Saturn are discussed in the Note after the Discussion section. A 
theoretical study of the distribution of objects in co-orbital motion is given 
by Morals & MorbideUi (2002). 

In this paper we describe a special set of collision orbits with the Earth 
which exist due to escape from L4 due to planetesimal perturbations. The 
perturbations cause a gradual peculiar velocity increase of the mass forming 
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at L4 so that it eventually achieves a critical escape velocity to send it to- 
ward a parabolic Earth collision approximately in the plane of motion of the 
Earth about the Sun. The region in velocity space where escape from L4 oc- 
curs in this fashion is relatively narrow. This mechanism therefore involves a 
special set of L4 ejection trajectories which creep towards collision with the 
Earth. At the end of Section 3 we will present a full simulation in three- 
dimensions of a collision of a Mars-sized impactor with the Earth assuming 
a thin planetesimal disk, using the general three-body problem, where plan- 
etesimal encounters with both the impactor and the Earth are done in a 
random fashion. The Appendix of this paper discusses the dynamics of the 
random planetesimal encounters. 

The paper has several main results: 
We show that a stable orbit at L4 exists where a Mars-sized giant impactor 
could grow by accretion without colliding with the Earth. We show that 
eventually perturbations by other planetesimals can cause the giant impactor 
to escape from L4 and send it onto a horseshoe orbit and then onto a creeping 
chaotic trajectory with an appreciable probability of having a near parabolic 
collision with the Earth. In the Discussion section we argue how this scenario 
fits in extremely well with giant impactor theory and explains the identical 
oxygen isotope abundances of the Earth and the Moon. The solar system 
itself provides a testing ground for our model. As we have mentioned the 
Trojan asteroids show that planetesimals can remain trapped at Lagrange 
points, and in the "Note added" we point out that the system of Saturn's 
moons provide examples where the phenomenon we are discussing can be 
observed, supporting our model. Finally both in the "Note added" and in 
the Appendix we discuss prospects for future work. 

The spirit of this paper is to suggest the intriguing possibility that the 
hypothesized Mars-sized impactor could have originated at L4(orL5). It is 
hoped that this lays the ground work for more detailed simulations and work 
in the future. 

2 Models and a Stability Theorem 

Let Pi represent the Sun, P2 the Earth, and P3 a third mass particle. We 
will model the motion of P3 with systems of differential equations for the 
restricted and general three-body problems. 

The first preliminary model, and key for this paper, is the planar circular 
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restricted three-body problem which assumes the following: 1. Pi, P2 move in 
mutual Keplerian circular orbits about their common center of mass which is 
placed at the origin of an inertial coordinate system X, Y. 2. The mass of P3 
is zero. Thus, P3 is gravitationally perturbed by Pi,P2, but not conversely. 
Letting rrik represent the masses of Pk, k = 1,2,3, then ms = 0, and we 
assume that 777,2/^1 = .000003. Let u be the constant frequency of circular 
motion of mi and 1712, uj = 2tt/P, where P is the period of the motion. We 
consider a rotating coordinate system (x, y) which rotates with the constant 
frequency lj as Pi and P2. In the x-y coordinate system the positions of Pi 
and P2 are fixed. Without loss of generality, we can set = 1 and place Pi 
at (/i, 0) and P2 at (— 1 + /i, 0). Here we normalize the mass of mi to 1 — /i 
and 7772 to /i, /X = 7772/ (7771 + 7772) = .000003 . The equations of motion for P3 
are 

X — 2y = X + fix 

(1) 

y + 2x = y + fly, 

where-^|,fi. = f, 

ri r2 

ri = distance of P3 to Pi = [{x — /i)^ + 7/^]^, and r2 = distance of P3 to 
P2 = [(x + 1 — /i)^ + y^] 2 , see Figure lU The right hand side of ((T)) represents 
the sum of the radially directed centrifugal force Fc = (x, y) and the sum 
Fg = {flx,f^y) of the gravitational forces due to Pi and P2. We note that 
the units of position, velocity and time are dimensionless. To obtain position 
in kilometers, the dimensionless position {x,y) is multiplied by 149,600,000 
which is the distance of the Earth to the Sun. To obtain the velocity in 
km/s, s = seconds, the velocity x,y is multiphed by the circular velocity of 
the Earth about the Sun, 29.78 km/s. For (Q), t = 2ti corresponds to 1 year. 

It is noted that is invariant under the transformation, x — x, 7/ — >■ 
—7/, t —t. This implies that solutions in the upper half-plane are symmet- 
ric to solutions in the lower-half plane with the direction of motion reversed. 
This implies, as noted in the introduction, that all the results we will obtain 
for L4 are automatically true for L5, and thus L4 need only be considered. 

System of differential equations has five equilibrium points at the well 
known Lagrange points L^, /c = 1, 2, 3, 4, 5, where x = 7/ = and x = 7/ = 0. 
(This implies that Fc+Fg = 0.) Placing P3 at any of these locations implies 
it will remain fixed at these positions for all time. The relative positions of 
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Figure 1: Rotating coordinate system and locations of the Lagrange points. 

Lk are shown in Figure H The locations of the Lagrange points for arbitrary 
G (0, 1) are a function of /i. Three of these points are collinear and lie 
on the and the two that lie off of the x-axis are called equilateral 

points. We note that the labeling of the locations of the Lagrange points 
varies throughout the literature. We are using the labeling consistent with 
that in Szebehely (1967), where in Figure is interior to to m2 and mi, 
and where L4 lies above the x-axis. (Note: This means L4 is 60° behind the 
Earth.) 

The three collinear Lagrange points Lkik = 1,2,3 lying on the x-axis 
are unstable. This implies that a gravitational perturbation of P3 at any of 
the collinear Lagrange points will cause P3 to move away from these points 
as time progresses since their solutions near any of these points are domi- 
nated by exponential terms with positive real eigenvalues (Conley 1969). The 
two equilateral Lagrange points are stable, so that if P3 were place at these 
points and gravitationally perturbed a small amount, it will remain in motion 
near these points for all time. This stability result for L4, L5 is subtle and 
was a motivation for the development of the so called Kolmogorov-Arnold- 
Moser(KAM) theorem on the stability of motion of quasi-periodic motion 
in general Hamiltonian systems of differential equations (Arnold 1961; 1989; 
Siegel & Moser 1971). A variation of this theorem was applied to the sta- 
bility problem of L^jL^ by Deprit & Deprit-Bartolome (1967). Their result 



7 



is summarized in the following result and represents a major application of 
KAM theory, 

L4, L5 are locally stable if < /i < /ii, fii = |(1 — |v^69) ~ .0385, and 
7^ = 2,3,4 fi2 = i(l - ^v^TSSS) ^ .0243, /xg = ^(1 - 1^^213) ^ 
.0135, /i4 ^ .0109. 

For further details connected with this result, see (Belbruno & Gott 2004) 
In our case, the Earth has /i = .000003 which is substantially less than fii 
and the exceptional values Hk, = 2, 3, 4 so that L4 is clearly stable for the 
case of the Earth, Sun system. 

An integral of motion for (Q) is the Jacobi energy given by 

J = -{x^ + f) + {x^ + y^) + /i(l - /i) + 2Vl. (2) 

Thus S(C) = &M.^\J = C, C G M} is a three-dimensional 

surface in the four-dimensional phase space (x, y, x, y), such that the solutions 
of (0) which start on S(C) remain on it for all time. C is called the Jacobi 
constant. The manifold S(C) exists in the four-dimensional phase space. 
It's topology changes as a function of the energy value C . This can be seen 
if we project S into the two-dimensional position space {x,y). This yields 
the Hill's regions 'H{C) where P3 is constrained to move. The qualitative 
appearance of the Hill regions 'H{C) for different values of C are described 
in Belbruno (2004) and Szebehely (1967). As C decreases in value, P3 has a 
higher velocity magnitude at a given point in the (x, |/)-plane. 

In this paper we will be considering cases where C is slightly less than 3, 
C ~ 3, where the Hill's region is then the entire plane. Thus, in this case P3 
is free to move throughout the entire plane. 

In the next section we will initially use the planar circular restricted three- 
body problem to obtain insight into the motion of P3 near L4. Ultimately 
we are interested in the general three-dimensional three-body problem with 
the mass points P^, k = 1, 2, 3 of respective masses m^. Unlike the restricted 
problem, m3 need not be zero, and Pi, P2 are not defined by constant circular 
Keplerian motion. Instead, Pi, P2 will be given initial conditions for uniform 
circular Keplerian motion between the Earth, P2, and Sun, Pi, assuming the 
Earth is lAU distant from the Sun. However, for 7^ this circular motion 
will not be constant. For 777.3 small, the deviation of the motion of Pi, P2 from 
the circular motion will in general be very small. Later in this paper we will 
be setting 

m3 = .l?7l2, (3) 



8 



where m2 is the mass of the Earth, and so /x = .000003. Thus, is a 
Mars-sized impactor. 

The differential equations for the general three-dimensional three-body 
problem in inertial coordinates (Xi, X2, X3) are defined by the motion of the 
3 mass particles Pk of masses rrik > 0, = 1, 2, 3, moving in three-dimensional 
space Xi, X2, under the classical Newtonian inverse square gravitational 
force law. We assume the Cartesian coordinates of the k-th particle are given 
by the real vector = {Xki, Xk2, Xks) G M^. The differential equations 
defining the motion of the particles are given by 

mfcXfc = > ^ ^ , (4) 

jit 



k = 1,2,3, where rj^ = \Xj — X^l = yYl^=ii-^ji ~ ^ki^ is the Euclidean 
distance between the fc-th and j-th particles, G is the universal gravitational 
constant, and ■ = J^. Equation (jU expresses the fact that the acceleration 
of the fc-th particle Pk is due to the sum of the forces of the 2 particles 
Pi, i = 1,2,3, i ^ k. The time variable t G M}. Without loss of generality, we 
place the center of mass of the three particles at the origin of the coordinate 
system. 

We note that the stability result by Deprit & Deprit-Bartolome (1967) 
provides conditions for stability of P3 with respect to L4; however, it does 
not necessarily provide conclusions on instability. It is proven more generally 
in Siegel & Moser (1971) that in the three-body problem, if satisfies 

27(mim2 -|- m2m^ + m^mi) > (mi + m2 + m^Y 

then the motion is unstable. This is not satisfied in our case since mi = 
1 — yU, 1712 = fi, ^3 = .Ifi, fi = .000003. However, if it is not satisfied 
it does not guarantee stability, and a deeper analysis is required such as 
KAM theory. In this more general result like that of Deprit & Deprit- 

Bartolome (1967) is not available. 

We have verified in the general three-dimensional three-body problem 
defined by (jH) with mi = 1 — /i, m2 = fi, = .l/x, /i = .000003, and 
more generally in the three-dimensional model for the solar system that L4 
is stable for a numerical integration time span of 10 million years. The model 
of the solar system we used includes the nine planets and is modeled as an 
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n-body problem with circular coplanar initial conditions using the current 
masses of the planets and radii of the planetary orbits. The integration time 
span of 10 million years is suitable for the purposes of our analysis. 

3 Chaotic Creeping Orbits Leading to Parabolic 
Earth Collision 

While we have shown via full solar system that L4 is stable it might be 
argued that this stability could be perturbed by other planetesimals, and in 
fact it is exactly this process that we are investigating (see Appendix). We 
expect that gravitational perturbations from other planetesimals will, via 
a random walk process in peculiar velocity, cause the Mars-sized impactor 
to eventually escape from L4. We numerically demonstrate in this section 
that there exists a family of trajectories leading from L4 to parabolic Earth 
collision. Producing these trajectories shows that Earth collision is likely 
when P3 escapes L4. P3 escapes L4 once it achieves a critical peculiar velocity 
- in the rotating frame. 

To describe the construction of the parabolic Earth colliding trajectories, 
we will begin first with the planar restricted problem. Then, we will show 
that the results hold up as we make the model more realistic. 

Assuming = we consider System ((T)) and place P3 precisely at L4. 
As long as the velocity of P3 relative to L4 is zero, then P3 will remain at L4 
for all time. 

The velocity vector at L4 for P3 is given by v = {x, y). Let a G [0, 27i] be 
the angle that v makes with the local axis through L4 that is parallel to the 
X-axis. Thus, v = V * (cos a, sin a), y = |v| = a/x^ + y'^. 

When V and if t = is the initial time for P3 at L4, then for t > 0, P3 
need not remain stationary at L4. If is sufficiently small, then by the result 
of Deprit & Deprit-Bartolome (1967) the velocity of P3 should remain small 
for alH > 0, and P3 should remain within a small bounded neighborhood of 
L4. This follows by continuity with respect to initial conditions. However, 
as ^(0) = V{t)\t=Q increases, then the resulting motion of P3 need not stay 
close to L4 for t > 0. This is investigated next. 

We fix a and fixing P3 at L4 at t = 0, we gradually increase ^(0) and 
observe the motion of the solution curve 7(t) = {x{t),y{t)) for t > for each 
choice of ^(0). This is done by numerical integration of System (^. (All the 
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numerical integrations in this paper are done using the numerical integrator 
NDSolve of Mathematica 4.2 until further notice.) The following general 
results are obtained which we first state, and then illustrate with a number 
of plots (In all of the plots of orbits of the restricted problem (P) in the x, y 
plane which are labeled 'Sun centered' the translation x—>a; + /i, has 
been applied which puts the Sun at the origin, and the Earth at the point 
(-1,0)). 

Rl For each choice of a G [0, 27r], as V^(0) is gradually increased from V{Qi) = 
0, and where 7(0) = (x(0),?/(0)) is at L4, the trajectory 7(t) for t > 0, 
remains in small arc-like regions about L4, which as ^(0) increases, evolve 
into thin horseshoe regions containing L4 and lying very near to the Earth's 
orbit about the Sun. As ^(0) increases further, the horseshoe region begins 
to close on itself, approaching forming a continuous annular ring about the 
Sun, coming close to connecting at the Earth. It is found that there exists 
a well defined critical value of V^(0) = V^*(0) where the ring closes at the 
Earth, and then the motion of P3 bifurcates from a motion constrained to 
the horseshoe-like region where it never makes a full cycle about the Sun, 
to a motion where it continuously cycles about the Sun, repeatably passing 
close to the Earth, and no longer in the horseshoe-like motion. ^*(0) has 
an approximate value for most values of a between .200 km/s to .600 km/s. 
We refer to this continuously cycling motion for V{0) = V*{0) as breakout. 
Breakout continues to occur for V{0) > V*{0) 

R2 In breakout motion for V{0) ~ V*{0) or 1/(0) = V*{0), the trajectory -f{t) 
for t > traces out a dense set of orbits in a thin annular region repeatably 
passing near the Earth, where the fiy-bys at Earth periapsis appear to be 
all approximately parabolic. The breakout orbit is chaotic in nature so that 
small changes in V{0) result in breakout trajectories which are in general 
significantly different in appearance, and still restricted to a thin annular 
region about Pi. The breakout trajectories as they cycle about the Sun 
have a high likelihood of colliding with the Earth. Moreover, for each a a 
near parabolic collision trajectory is readily found for V{0) ~ V*{0). The 
collision orbits move near the Earth's orbit and gradually approach the Earth 
for collision. (This gradual motion approximately along the Earth's orbit, 
we refer to as creeping.) The collision orbits can creep to collision along the 
direct or retrograde directions with respect to Pi. 

Demonstration of Rl 
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We choose an arbitrary velocity direction v for P3 at L4 at t = 0, where 
V points in the vertical positive y direction. In our simulations, the location 
of L4 is (—.5, and the y— coordinate is input with the value .866025404. 
Beginning with a = 7i/2, we choose a magnitude V = V{0) = .001, and 
numerically integrate the system of differential equations forward for 
t G [0,1000]. This velocity magnitude is small, and since L4 is stable P3 
remains in a thin arc-like region approximately of radius 1 shown in Figure |21 
P3 starts at the location x(0) = (.5, .5v^), and moves down in the posigrade 




Figure 2: 'y(t),V{0) = .001, t G [0,1000], x vs y (i.e. x-axis is horizontal, 
y-axis is vertical). Sun centered. 

direction with respect to the Sun. As it moves, it performs many small loops 
as are shown in Figure |21 These loops occur since the semi-major axis of 
the orbit of P3 has changed slightly from 1 and the orbit of P3 has a slight 
nonzero eUipticity, both due to the addition of V{0). So, as it moves in 
its approximate elliptical motion over the course of one year it falls slightly 
behind and forward with respect to the Earth when it is at its apoapsis and 
periapsis, respectively. Each loop forms in one year. Thus, for t = 1000, 
there are 1000/(27r) loops. P3 moves down to a minimal location where 
y is approximately .7, and then it turns around and moves in the upward 
direction where the small loops point in the opposite direction when it was 
moving in the downward direction. The superposition of the loops makes 
a braided pattern as seen in the lower half of Figure |21 P3 stays in this 
bounded arc-like region since L4 is stable, and the velocity V{0) is relatively 
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small. (If V{0) = 0, then P3 stays fixed at L4 for all time.) Because the 
velocity magnitude is small, P4 has a Kepler energy nearly that of L4, and so 
its semi-major axis with respect to the Sun deviates from 1 by a negligible 
amount. Thus, as it moves, it stays nearly on a circle of radius 1. That 
is, in an inertial coordinate system it stays approximately on Earth's orbit 
about the Sun. As long as V{0) is small, which it is throughout this paper, 
the trajectories of P3 remain close to the Earth's orbit and move with small 
loops, in the rotating coordinate system. The particle P3 creeps slowly along 
the Earth's orbit initially in a posigrade fashion, and then in a retrograde 
fashion away from the Earth. 

The above procedure is repeated, where we slightly increase the value of 
V^(0) to .004 at L4 at t = 0. Since V^(0) has increased, then as is seen in Figure 
EJ where t G [0, 1000], P3 creeps further in its Earth-like orbit about the Sun. 




Figure 3: 7(t), 1/(0) = .004, t e [0, 1000], x vs y. Sun centered. 

Since V{0) is small, the trajectory of P3 deviates slightly from a circle of 
radius 1. This deviation slightly increases as V{0) increases. The addition of 
V{0) at L4 causes P3 to have a slightly smaller value of the Jacobi integral, to 
be slightly less than 3 (C ~ 3). This means that P3 becomes more energetic, 
and thus can creep further along the Earth's orbit. Increasing V{0) by .001 
to .005 causes the increased creeping shown in Figure El where t e [0, 1000]. 
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Figure 4: 7(t), V{0) = .005, t e [0, 1000], x vs y. Sun centered. 




Figure 5: 7(t), 1/(0) = .009, t G [0, 1000], x vs y. Sun centered. 

In Figure 1^(0) is increased to .009. P3 leaves L4, moves downward 
in a posigrade fashion to slightly behind the Earth, then turns around and 
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moves in a retrograde fashion on its Earth-like orbit about the Sun, until it 
approaches the Earth from the front turning around and then moving in a 
posigrade fashion. 

A braided pattern results due to the fact the Earth-like orbit is traversed 
twice, with loops pointing in the inner and outer directions. The resulting 
complicated looking trajectory is symmetric with respect to the x-axis due 
the symmetry mentioned earlier for the restricted problem. The width of the 
region near the Earth's orbit in which P3 moves has slightly increased due to 
the increase in V{0). 

We note that general appearance of the orbit of the asteroid 2002 AA29, 
mentioned in the introduction (Conners et al. 2002), remarkably is very sim- 
ilar in appearance to Figure El Unlike the planar orbit considered here, 2002 
^4^429 has a inclination of 10 degrees with respect to the plane of the Earth's 
orbit. It's oscillation period is approximately 95 years. The approximate 
period of the orbit in Figure El is about 159 years, which is not too dissimi- 
lar. Orbits of this type are called horseshoe orbits. The horseshoe orbits are 
constrained to a region we refer to as a horseshoe region so P3 cannot move 
past the Earth. Such a region is discussed in Murray & Dermott (1999). Let 
6 be the polar angle measured from the positive x-axis for the position of 
P3. The horseshoe orbits have the the property that 6 ^ it. This means that 
P3 will not fly by the Earth. For other papers on this motion, see (Christou 
2000; Hollabaugh & Everhart 1973; Mikkola & Innanen 1990; Namouni 1999; 
Weissman & Wetherill 1974). 

When ^(0) reaches .011, P3 is able to escape from the thin horseshoe-like 
region and fly by the Earth as is seen in Figure IHl 

This achieves breakout motion where P3 then cycles about the Sun only in one 
direction. In Figure El the cycling is in the retrograde direction. This actual 
cycling is not shown in this figure since for the time range given, breakout 
into cycling motion occurs when t = 988, on the outer retrograde trajectory. 
P3 first leaves L4 moves near the Earth, then back up in a retrograde fashion 
going all the way around the Sun to near and in front of the Earth, then 
moving around the Sun again in a posigrade fashion to its location behind the 
Earth, then finally it moves on the outer trajectory in a retrograde fashion 
back to just ahead of the Earth when it crosses by the Earth at t = 988 
(crossing the x-axis near the Earth), then performing the cycling breakout 
motion after that time. This transition from creeping horseshoe motion to 
creeping breakout motion is what is desired for this paper. The transition 
from horseshoe motion to breakout motion represents a bifurcation from one 
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Figure 6: j{t), V{0) = .011, t e [0, 1000], x vs y. Sun centered. 



type of motion to a different type. We are interested in tiie likelihood of 
Earth collision while in breakout motion just after the bifurcation. This 
represents breakout motion with minimal energy. 

This transitional breakout motion has two important properties: 

1. P3 moves in a thin annular region about the Sun, 

2. P3 repeatably flys by the Earth. 

These properties imply the following: Since the annular region is thin, the 
Earth fly-bys are in general close. The close Earth fly-bys are approximately 
parabolic in nature, as we will demonstrate, and as P3 flys by the Earth its 
actual velocity vector is approximately tangent with the Earth's orbit. This 
implies that P3 gains a negligible velocity increase due to gravity assist as 
it flys by the Earth as we will show. This guarantees that P3 will continue 
to move in an Earth-like orbit about the Sun, and continue to cycle. This 
implies that as P3 moves around the Sun, it will densely flll the thin annular 
region it moves in. This means that it has a high likelihood of colliding with 
the Earth. We will demonstrate that coUision readily occurs in these creeping 
breakout orbits. 

In fact, the previous case where V{0) = .011, which is the flrst breakout 
motion we computed, leads immediately to collision at t = 1384.7176 (or 
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220.3847 years). In our exposition below, we will use a slightly different 
value of V^(0) which happens to achieve collision at an even earlier time. 

Breakout motion is seen in Figure [7| It is observed that a shift from 
V{0) = .009 to ^(0) = .012 causes a qualitatively different looking picture, 
where the bifurcation between horseshoe and breakout motion is clearly seen. 




Figure 7: 7(t), 1/(0) = .012, t E [0, 1000], x vs y. Sun centered. 

The case just considered is for the direction a = tt/2. The same procedure 
produces critical values of V{0) = V{0)* leading to breakout motion, from 
horseshoe motion, for any value of a; G [0, 2tt]. A set of these for a increments 
of 7r/8 are listed in (Belbruno & Gott 2004) in Table 1. This is graphically 
shown in Figure |H1 In this figure, the length of each line is equal to the value 
of V^*(0) in that direction. In this way, a smooth variation of V*{0) as a 
function of a is numerically obtained. There is a sharp spike in the value 
of V*{0) which has a maximum at 5.1027r/8 of .22. There is also a similar 
maximum near the value of 13.57r/8. These are not listed since they are not 
typical: almost all the values of V*{0) are in the range of values illustrated. 
The minimum value of V* (0) = .0057 is for a = (97r/8) — .01. Multiplying the 
values of l^*(0) by 29.78 yields a range of velocity values generally between 
180 meters/s and 1.2 km/s. Note that the two directions corresponding to 
the maximal spikes in velocity seen in Figure |H1 approximately lie near the 
Sun and anti-Sun directions. In this figure the Sun is toward the lower right. 
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Figure 8: Initial velocity directions v(0) at L4 whose magnitude corresponds 
to the associated critical breakout velocity V*{0), which are hsted in Table 
??. 

Note that the method, or algorithm, used above to estimate the critical 
velocities V*{0) at L4 leading to breakout motion, is similar in nature to 
the method of estimating transitional stability regions, called weak stability 
boundaries, between capture and escape about the Moon described in Bel- 
bruno (2004). This capture region has important applications. It was used 
by one of us to find a new type of low energy route to the Moon in 1990 
where lunar capture is automatic (Belbruno & Miller 1990). This special lu- 
nar transfer was designed in order to resurrect a Japanese lunar mission and 
enable the spacecraft Hiten to successfully reach the Moon in October 1991 
with almost no fuel (Belbruno 1992; Frank 1994). More general references 
on this are (Adler 2000; Belbruno 2004; Belbruno & Miller 1993). 

Other methods could be used to study the bifurcation from horseshoe to 
breakout motion such as the computation of suitable surfaces of section to the 
trajectories in phase space, and then monitoring the iterates of intersecting 
trajectories on the section. This would give a more complete knowledge of 
the phase space near breakout motion, but this approach is not necessary for 
our purposes. The algorithm we have described accurately determines when 
bifurcation occurs. 

We have also performed an analysis to understand the relationship of the 
critical breakout velocities as a function of the mass of the Earth m2. It is 
found that roughly 

V*{0) oc ml^^. 

As an example, we consider the two cases: m2 = .Im^, .Olm^, = 
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.000003. When ma = .Ittie, then for a = 0, 7r/2, vr, 37r/2, we obtain V*{0) = 
.004, .007, .004, .007, respectively, and for = .Olm^, we obtain V*{0) = 
.002., 003, .002, .003, respectively. These results imply that V*{0) ~ 0.6ml^^ 
for a = 0,71, and when a = 7r/2,37r/2, V^*(0) ~ ""^a^^- (Thus a giant 
impactor trapped in a stable orbit about L4 and unperturbed will remain 
trapped there as the proto-Earth grows by accretion. Breakout velocity in- 
creases as the proto-Earth grows, postponing breakout, but at late times 
after the proto-Earth has reached essentially its full mass, according to our 
scenario, perturbations can drive the giant impactor to breakout.) 
This concludes the demonstration of Rl. 

Demonstration of R2 

We first show how to readily find trajectories from L4 which collide with 
the Earth. The value of a = 7r/2 is again considered, and we consider the case 
y(0) = .012 ~ V*{^) shown in Figure [7| Plotting the distance r2 between 
P3 and the Earth for t G [0, 1000] reveals the times of the various Earth 
fly-bys. It was found that the case of ^(0) = .012, for the given range of 
had very close Earth fiy-bys, but no actual collision. Randomly altering this 
value of V^(0) yielded a collision on our second random choice of values of 
Vifi) = .0119981. This is seen by plotting r2 as a function of time shown in 
Figure M 

By magnifying the regions near minima of r2 it can be seen which ones 
may yield collision. In this case, counting from the left to the right, we 
determined that the first, fourth, fifth minima yield distant fiy-bys at over 1 
million km. Thus, in these cases, the Earth fiy-bys are not of interest. The 
third fiy-by misses the Earth by about 18,000 km. However, it is found that 
the second fly-by in fact collides with the Earth. The time of collision with 
the surface of the Earth is at t = 360.181558, corresponding to 57.3247 years. 
This time is calculated when the center of the impactor, viewed as a circle of 
radius rj = 3397 km, intersects the surface boundary of the Earth, which is 
at a radial distance te = 6378.14 km, from the Earth's center. This is seen 
in Figure ITUl 

The horizontal line indicates the Earth's radius of 6378.14 km here repre- 
sented by the value r2 = .0000426346. The time of collision with the surface 
of the Earth is at t = 360.181558, corresponding to 57.3247 years. We are 
assuming that each point of the trajectory of the impactor at any given time 
is at the center of the impactor. More accurately, however, collision actu- 
ally occurs a few moments earlier when the surface boundary, circle, of the 
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Figure 9: Variation of distance r2 of P3 to P2 (earth) as a function of i e 
[0, 1000] in dimensionless units, ^(0) = .0119981. 
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Figure 10: Distance r2 of P3 to P2 (earth) goes below the earth's radius 
.0000426346(i.e. 6378.14 km) proving coUision has occurred. V{0) — 
.0119981. Plot is t vs r2 in dimensionless units. 

impactor touches the surface boundary, circle, of the Earth. That is, when 
r2 = rj + te = 9775.14 km, where for rj we take the radius of Mars, since 
this is a Mars-sized object, or, in dimensionless units, when r2 = .00065342. 
In general, we assume collisions mathematically occur in our numerical sim- 
ulations when r2 < .00065342. 
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We now show what the colhsion trajectory looks hke and discuss its prop- 
erties. The colhsion trajectory, CI, is shown in Figure ITTl 




Figure 11: Entire collision orbit (x vs y)- Originating at L4 at t = and 
colliding with the earth when t = 360.18 or, equivalently, 57.32 years, V{0) = 
.0119981, Sun centered. 

It starts at L4, moves in a posigrade fashion toward the Earth, turns 
around and in a retrograde motion moves around the Sun to collide with the 
Earth. A view of this orbit in its final 9.57 years is shown in Figure IT^ 

Collision with the Earth itself and the final 9.14 hours of the trajectory 
are shown in Figure IT^ 

In this figure we stopped the trajectory of the impactor before its Earth 
periapsis. However, if it were continued beyond collision it would reach 
its periapsis point approximately on the x-axis inside the Earth's radius at 
t = 360.1817 where rs = .00003. (See Figured). 

If it were extended so that t G [0, 1000], the trajectory would be as shown in 
Figure IT^ 

This figure for V{0) = .0119981 is shown to compare with Figure[7|for V{0) = 
.012 indicating the sensitive, or chaotic, nature of the breakout motion, where 
a difference of V^(0) by .000002 yields a qualitatively different appearing 
trajectory. The chaotic nature of the motion near breakout is also seen in 
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Figure 12: Orbit approaching collision with the earth (x vs y). Time duration 
of 9.57 years shown, t G [300,360.181558], axis earth centered. 

Figure El which is a plot of r2 for V^(0) = .0119986, when compared to 
FigureHwhere 1/(0) = .0119981. 

It is seen that although the difference in V{0) = .0000005 at L4, there is 
a significant difference in the qualitative appearance of the two plots. This 
is caused by the fact that infinitesimally small changes in V{0) at L4 can 
cause slightly different Earth fly-by conditions which over long time spans, 
t G [0, 1000] can cause the trajectory to change noticeably if any of the fly- 
bys are close. However, as we will see in the following, close fiy-bys will 
only yield negligible Kepler energy increases with respect to the Sun. So, 
although the trajectory may have a qualitative different appearance, it will 
still have approximately the same Kepler energy before and after close fiy-bys. 
The change in the trajectories for tiny changes in V{0) observed is typical 
for chaotic motion in general, and is a sign that a hyperbolic invariant set 
likely exists in the phase space for the breakout motion of P3. A hyperbolic 
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Figure 13: Collision of the impactor with the earth, x vs y. Final 9.14 hours 
of trajectory shown, earth centered. 

invariant set in general is a Cantor set all of whose points are have a saddle- 
like structure produced by a transverse homochnic orbit, whose existence is 
given by the Smale-Birkhoff theorem (Bclbruno 2004). The use of the term 
chaotic in a strict mathematical sense means the existence of a hyperbolic 
invariant set. 

It is remarked that CI represents a physical collision with the surface 
of the Earth, where, at periapsis below the Earth's surface, r2 = .00003. 
It turns out that actual pure collisions where r2 ~ to high precision are 
readily found as well. For example, V{0) = .011998 leads to a pure collision 
at t — 360.1898. In this type of collision, P3 asymptotically approaches the 
collision manifold which is a set of measure zero. 

It is observed that since the motion of P3 repeatably passes near to the 
Earth in breakout motion, the Earth tends to readily pull P3 toward pure and 
physical collisions. The set, or manifold, of pure collision trajectories are a 
subset of physical collision trajectories, and, in fact, are a set of measure zero 
in the four-dimensional phase space of position and velocity (Belbruno 2004). 
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Figure 14: Continuation of collision orbit across the a;- axis when t = 
360.1816, X vs y, earth centered. 




Figure 15: Extended collision orbit, through collision, t e [0, 1000] ( i.e. 
t e [0, 159.15] years), x vs y. Sun centered. 
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Figure 16: Variation of r2,t e [0,1000], for 1^(0) = .0119986. Compare with 
Figure 13 where 1^(0) differs by .0000005. 

Since they are a set of measure zero, their near occurrence is reflective of the 
fact the fly-bys of the Earth are close and that the Earth has a considerable 
gravitational focusing effect when the trajectory is near parabolic. 

It turns out that CI is approximately parabolic at collision. This is seen 
by plotting the Kepler energy E2 of CI with respect to the Earth. In inertial 
Earth centered coordinates, X = {Xi,X2), E2 = ||Xp — j^. In barycentric 
rotating coordinates x = {xi,X2) = {x,y), E2 is transformed into 

E.^l\M'--^^^lrl-L, (5) 

where, r2 = + 1 — /i)^ + L = X1X2 — X2(xi + 1 — /i), (Belbruno 

2004, 2002). L is the angular momentum of P3. ®, is evaluated along 
CI and plotted in Figure IT7I as a function of t G [0, 360.181558]. From Figure 
ITS| E2 = .000054 at collision, which is nearly parabolic. 
(If the collision were purely parabolic then E2 = 0.) This yields a very slight 
hyperbolicity whose hyperbolic excess velocity with respect to the Earth, 
Voo = V2E2, has the value of .0104. This is equal to V*{0) for a = 7r/2 to 
within .0014. In scaled coordinates, V^o = 310m/s and V*{0) = 328m/s. 
Voo is close to V*{0) because at Earth periapsis V*{0) is like the velocity 
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Figure 17: Kepler energy E2 with respect to P2 as a function of t along entire 
collision orbit. It is seen that £'2 — > 0. 

at infinity. This velocity is approximately maintained along the orbit as it 
approaches collision. At actual Earth fly-by at periapsis the velocity with 
respect to the Earth increases due to the attraction of the Earth, and for this 
collision orbit V = .44 at Earth periapsis. 

In (Belbruno & Gott 2004) it is analytically shown that in the critical 
or near critical breakout motion, all close Earth fiy-bys, including collision 
trajectories, are approximately parabolic at periapsis. For critical breakout 
trajectories, which start at L4 at time t = 0, V{0) = V*{0). For near critical 

breakout motion we assume that V{0) - ^^^(O). Notationally, V{0) - V*{0) 
includes both of these cases. 

We define the terms, 'close Earth fly-by', 'approximately parabolic'. Let 
j{t) be a trajectory which performs a fly-by of the Earth, with a periapsis 
distance r2 at some time t. We say that this is a close Earth fly-by if r2 < 
100, 000 km, or in dimensionless coordinates, r2 < .000668. The figure of 
100,000 km is arbitrarily chosen since for weakly hyperbolic fiy-bys of the 
Earth beyond this distance, the effect of an Earth gravity assist is neghgible. 
Physical collisions are included as close fiy-bys. 

We use E2 to determine the type of coUision, which is computed at Earth 
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Figure 18: E2 as a function of t at the end of the colhsion orbit approaching 
the value of .0000541 at periapsis below the earth's surface, where r2 = .00003 
(i.e. 4488 km), for t — .00014 (i.e. 12 minutes) beyond earth coUision, 
t e [360.175,360.1817]. 

periapsis. So, in the case of physical collision at the Earth's surface, we 
propagate the trajectory to Earth periapsis within the Earth. This point 
occurs at a very short time after physical collision which for the case of C/ is 
only 12 minutes. At the periapsis point E2 is evaluated at the trajectory state 
of position and velocity. If |ii^2| ~ 0, the collision or coUision trajectory is 
called approximately parabolic. It could be slightly elliptic, slightly hyperbolic 
or purely parabolic. It turns out, as we will see, that for break out or near 
breakout motion, the fiy-bys are all approximately parabolic. The following 
result is obtained (For details see Belbruno & Gott 2004): 

For the set of critical or near critical breakout velocities at L4, the value of 
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i?2 at the close Earth fly-bys at periapsis has the value 

E,^^ViOf^^V*{Of^O; (6) 

That is, the close Earth fly-bys are approximately parabolic. This is true for 
all the values of a except those values in small neighborhoods o/5.1027r/8, 13.57r/8 
(see comment below). 

This implies that a trajectory 7(t) starting near critical breakout velocity 
at L4 for t = to will satisfy equationlHlfor any future time t > to corresponding 
to any close Earth fiyby at periapsis. 

More precisely, as shown in Belbruno & Gott 2004, for a trajectory start- 
ing at L4, at Earth periapsis on a close Earth fly-by at a distance r2 = 62, 

E2 = \v\Q) + Oo(/i) + + 02(^2) (7) 

where Oo = 3fi, d = 26i, O2 = 252Cos^2, r2 = 62 < .000668 < 1, n = 
1 + 5i, \Si\ < .000668, \5i\ > 0, ^2 > 0. This relation yields (0) for fx,5i,5i 
small. 

Let P3 be at L4(or L5) at t = 0, and let v(0) be the initial velocity 
with magnitude ^(0) = |v(0)|. Then, the Jacobi integral J has the value 
J = Co = 3 — V^^(O) — — jj). This implies that for the set of critical 
breakout velocities V*(0) at L4 for a G [0, 27r] (see Figure |HI), Co ~ 3. 

As mentioned earlier, there are two sharp spikes in the breakout velocities 
shown in Figure |H1 of values .22, .25 which occur for a = 5.1027r/8, 13.57r/8, 
respectively. However, most values of V*{0) vary between approximately .006 
and .05 if two intervals in a of total width approximately .157 radians are 
deleted near where the spikes occur. This implies that for nearly all of the 
values of a, equation ^ implies that approximately 

^2 e [.000018, .0012]. (8) 

Thus, P3 is approximately parabolic at collision. The range given by (jSI) 
is a crude estimate of E2. The observed value for CI of E2 = .000054 is 
contained within this interval. A sharper estimate can be made using ([7j). 
For CI, V{0) = .0119981, and at fly-by periapsis, below the Earth's surface, 
r2 = 62 = .00003. This occurs approximately on the xi— axis implying 
n = -61 = -.00003, 02 = 0. Substitution of V{0),ri,r2,92 into (0) yields 
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the value, E2 ~ .00007. Noting that the numerically observed value of E2 at 
periapsis for CI E2 = .000054, the predicted value is in error by only .000016 
demonstrating the accuracy of the predicted values. 

We remark that in all the cases numerically observed, E2 was not neg- 
ative(i.e. the orbit was not elliptic) during the fly-by. Some such fly-bys 
are likely to exist due to the chaotic nature of breakout motion, however, 
the probability of finding trajectories with elliptic fly-by states is apparently 
small. Their chaotic structure and low probability of occurrence is stud- 
ied in Belbruno (2004). When E2 is near to zero this defines weak capture 
studied in Belbruno (2004), where P3 will, in general, move about the Earth 
in a chaotic fashion generally leading to escape or collision. However, as 
remarked, the case of interest here is when E2 is very slightly hyperbolic. 

The retrograde collision trajectory CI emanating from L4 is paired with 
another symmetric collision trajectory CI* emanating from L5 which is sym- 
metrical to CI and moves in a posigrade fashion about the Sun. This follows 
by the symmetry of solutions mentioned in Section |21 It will collide with the 
Earth in the 4th quadrant as shown in Figure IT^ 

Probability of Collision at Breakout for the Restricted Problem 

A measure of the likelihood of finding collision trajectories is now de- 
scribed. This is done for the four basic initial velocity directions at L4: 
a = 0, 7r/2, TT, 37r/2. For an initial velocity of P3 at L4 for a given a we as- 
sume the corresponding breakout velocity V^*(0) as shown in Figure |H1 The 
orbit of P3 is propagated from L4 for t > and since it is in breakout motion 
we know that it will not be in horseshoe motion, but will cycle about the Sun 
and repeatably fly past the Earth. We can numerically demonstrate that col- 
lision with the Earth is likely. This intuitively makes sense since the fly-bys 
will be close and the the annular region supporting the breakout motion is 
narrow. Now, for a given initial velocity at L4 for t = we see from Figure |H1 
that V*{0) is given up to three digits. For a given value of V*{0), depending 
on a, we propagate the trajectory for up to t = 4000, which corresponds to 
637 years, and see if collision has occurred, t = 4000 is chosen arbitrarily, 
for convenience and is fairly small in astronomical terms. If no collision oc- 
curred in that time, then we give V*{0) a random perturbation by adding 
to it the random number .OOOmn, where m, n are positive random integers 
ranging from to 9. For a choice of m, n the trajectory is propagated again. 
If collision does not occur, we repeat the process again for a different choice 
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of m, n, continuing trials until success is achieved. 

For a = 0,7r/2, we required two random trials for success, where success 
means we achieve collision within t = 4000 = 637years. For a = n, three 
random trials were required until we have success, and for a = 37r/2, six 
random trials were required for success. Therefore we have achieved success 
in four random trials out of thirteen. This gives our best estimate of the 
probability V of success for < t < 4000 as 

P ~ 1. (9) 

If we had not limited ourselves to t = 4000 the probability would have been 
larger. This probability is discussed in further detail in the Appendix. We 
have run a sufficient number of trials to produce a rough order of magnitude 
estimate of this probability which is sufficient for our purposes, but a large 
number of additional trials could establish this number to higher accuracy. 

It is noted that the gravitational focusing on P3 to cause a collision is 
substantial. This is related to the fact the breakout motion is occurring at 
a fixed energy for the planar restricted problem. The fixed energy yields a 
three-dimensional energy surface obtained from the Jacobi integral. As is 
proven in Belbruno (2004), the manifolds leading to collision at P2 are two- 
dimensional, and although they arc a set of measure zero, the particle P3 is 
readily able to move asymptotically close to these surfaces and to collision 
after the gravitational focusing. The collision manifolds on the Jacobi integral 
surface separate the phase space, so it is fairly easy for P3 to get near to the 
collision manifold. In higher dimensions this separation of the phase space on 
the Jacobi surface does not occur, and the collision manifold is more elusive. 

This concludes the demonstration of R2. 

It is interesting that these creeping chaotic orbits seem to lead naturally 
to collision with Earth (as proposed by the giant impactor theory) rather 
than to capture into a bound orbit (as in the sister planet theory). If one 
wanted to have a sister planet theory, of course, L4 would be a promising 
place for the Moon to start out. So it is significant that our chaotic creeping 
orbits (shghtly hyperbohc-nearly parabohc) lead naturally to coUision rather 
than capture. This favors the giant impactor theory. The sister planet theory 
would of course also have a problem with the difference in iron between the 
Earth and the Moon. 

Random Walk, AV Accumulation, and Relevant V{0) Range 
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In determining V*{0) at L4 above, we kept P3 fixed at L4 and gradually 
increased V^(0) for a given velocity direction. This yields a well defined set 
of V*{0){a) for a G [0,2%]. 

We now consider a more realistic way that P3 would increase its veloc- 
ity in a gradual fashion. The mechanism for this is to assume that P3 is 
randomly being perturbed by encountering other planetesimals (whether by 
gravitational encounter or direct collision) and in each encounter, it acquires 
an instantaneous kick AV. So, it is not kept fixed at L4. To make this 
more realistic, we assume that the times of encounters are random, within a 
large range, and the direction a of the kicks are random. The only thing we 
normalize is the magnitude of the |Ay|'s which for convenience is held fixed. 

Thus, P3 starts at L4 with a zero velocity, and at time t = ti = a 
velocity V^(0) = AV^, is applied in a random direction. This yields a vector 
vi with magnitude AV. P3 moves on a trajectory 'y{t) in a neighborhood of 
L4, assuming that the value of AV is small. At a random time ^2 > another 
velocity vector V2 of random direction and magnitude AV is vectorially added 
to P3's velocity ait = t2. Then the trajectory is propagated for t > t2 until at 
another random time t^> t2 a random vector V2 of magnitude is vectorially 
added to Ps's velocity vector at t = t^, and this process continues creating a 
sequence tk of times, tk+i > tk, and velocities Vfe, A; = 1, 2, 3, . . . . 

While the AVs are being applied, the trajectory ^(t) is gradually moving 
further from L4, but since the velocity directions Vfc are applied randomly, 
the path of the trajectory 7(t) will move further away from L4 for some time 
spans, and then move toward L4 for others. However, as k increases, one 
would expect, by the principle of random walk, for P3 to eventually escape 
L4 and creep toward the Earth for k sufficiently large when the velocities 
Vfc, A; = 1, 2, 3, . . . applied on the trajectory ^{tk), gradually accumulate to a 
sufficiently large magnitude for breakout to occur. If the were all applied 
in the direction of motion of P3 at tk, then the magnitudes AV would add 
producing an cumulative velocity addition of kAV at the A;th step. However, 
the directions of Vjt are random, and by the principle of a random walk, the 
number of encounters k before ejection occurs should be expected to instead 
satisfy 

VkAV ^ .006 (10) 

for k sufficiently large(and AV sufficiently small) where .006 is approximately 
the minimum value .0057 of {V^*(0)}. This makes dynamical sense, since as 
the Ay's are applied, the trajectory 7(t) would seek to minimize the Jacobi 
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energy, and hence the velocity, along its path. We found in all our numerical 
simulations the following result: 

For a given value of AV , the number k of random applications required 
for breakout to occur approximately satisfies ilU^) . 

We describe this process of random walk AV accumulation, and its veri- 
fication. 

For convenience we choose AV^ = .001 and start at L4 with zero initial 
velocity. (fTUj) implies that k should satisfy \/kAV ~ .006, which yields 
k = 36. It was found that breakout occurred when A; = 36 as predicted. 
Random kicks of velocity AV were applied in random directions and after 
random time intervals t^- (As in all future runs the times between kicks are 
just chosen to be large enough to randomize the position. The real time for 
random walkout is expected to be much longer in years - perhaps 30 million 
years as considered in the Appendix.) 

From the above, we have the following result, 

(Random Walk AV^ Accumulation) Under a realistic assumption of random 
walk, the peculiar velocity for P3 accumulates proportional to the square root 
of the number of encounters until it reaches a breakout state. Since a random 
walk is isotropic the peculiar velocity is likely to encounter the breakout state 
first at a point near the minimum value of .006 of the set {^*(0)} thus giving 

Therefore, substituting V*{0) = .006 into © implies that for close Earth 
fiy-bys resulting from the random walk process at or near breakout, E2 ~ 
.000018. This implies that at close Earth fiy-by resulting from the random 
walk process, a nominal value of V^o = -006 which is 179 m/s. When P3 does 
a close fiy-by of the Earth, after passage through periapsis it will receive 
a gravity assist and increase, or decrease, its velocity with respect to the 
Sun. A measure of this velocity change is observed due to the bending 
of the trajectory of P3 as it passes through periapsis. For example, this 
bending is clearly seen in Figure The more distant the fiy-by, then in 
general the less the bending. The maximum bending is obtained from pure 
collision trajectories, where the bending angle is ivr. It is determined in 
Broucke (1994) that the resulting change in magnitude of the velocity, 6v 
with respect to the Sun due to gravity assist is maximally, 2Voo. Thus, for 
each close Earth fly-by, the expected maximum gain in velocity magnitude is 
approximately 358 m/s. In general, they will be less. 
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The maximal velocity of 358 m/s is a relatively small number and will 
have little effect on a breakout trajectory when it has a close Earth fly- by. 
This velocity is less than .012 of the orbital velocity, inducing eccentricities 
into the trajectory of P3 after fiy-bys of at most this order. It is found in 
general that within time spans on the order of 2000 time units, there are 
generally only one or two close Earth fiy-bys. This implies that P3 will 
remain in breakout motion about the Sun in a relatively thin annular region 
for very long periods of time, generally tens of thousands of time units, and 
repeatably pass by the Earth without being ejected. 
Collisions in Three Dimensions and the Mars Sized Impactor 

Thus far we have constrained P3 to lie in the plane of motion of the Earth 
and Sun in the planar restricted three-body problem. It that situation it was 
seen that collision trajectories are readily found. When the motion of P3 has 
an out of plane component, z, added to it, it is more complicated. In this 
case we have the three-dimensional circular restricted three-body problem, 
defined exactly the same way as for the planar problem, except that P3 can 
move in the z direction, see (Szebehely 69; Belbruno 2004). By continuity 
with respect to initial conditions, it must be the case that collisions in the 
planar case will persist in the three-dimensional case if |i;| is sufficiently small. 
A way to approximately get a measure of the maximal allowed z motion is 
to consider a collision trajectory from L4 to the Earth generated at or near 
critical breakout, and then see how much i can be added at L4 and still 
maintain collision with the Earth which is now a three-dimensional sphere. 
This is a relatively straight forward calculation. 

In the Appendix we show that Earth collisions should persist providing 
|i;| ~ .0034. This implies a thin disk of planetesimals. In the Appendix we 
show that a disk of thickness 28rE easily satisfies this requirement. This 
has an angular width of 4.1 minutes of arc as seen from the Sun. Although 
this seems thin, it turns out that inner B, A rings of Saturn, extending from 
92,000 km to 140,210 km, have a thickness of .18 - 1.7 seconds of arc as 
seen from the center of Saturn, which is even thinner than required in our 
situation. 

The modeling of most interest in this paper is for the general three- 
dimensional three-body problem defined by equation 0J where 7713 = .lm2 7^ 
0, so a Mars-sized Earth impactor is modeled. This is used by Canup & 
Asphaug (2001). Although the motion of the Earth, P2, is given initial 
conditions for uniform circular motion, about the Sun, Pi, it need not remain 
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circular as time progresses due to the gravitational perturbations of P3. This 
property makes the problem more interesting. To better understand this 
and to see its effect on collision trajectories, the general planar three-body 
problem is first considered. It is defined from Q by setting X^s = 0,k = 
1,2,3. 

As with the restricted problem, a rotating coordinate system {x, y) is 
chosen, this time, rotating with the mean motion of the Earth about the 
Sun. In this system the Earth is not fixed on the -x axis as in the restricted 
problem since it is perturbed by P3 especially during fly-bys. 

Just as in the restricted three-body problem, breakout from L4 can be 
defined for the general planar three-body problem with exactly the same 
methodology as for the restricted problem, by gradually increasing the ve- 
locity of P3 at L4 until breakout is achieved, for any given velocity direction. 
Nearly identical results are obtained as in Rl (Figure |H1) and this analysis is 
not duplicated for this situation. The desired random walk AV accumula- 
tion process is defined exactly as before for this current problem, and we 
have verified that the same results are obtained. 

We determine the breakout of P3 to Earth collision as we did in the 
restricted problem by gradually increasing the velocity at L4 until a critical 
velocity is reached. That is, we are not modeling the random walk process, 
and are copying the procedure we performed for the restricted problem. This 
is done to keep the modeling as close as possible to the restricted problem 
in order to better understand how P3 perturbs P2 and the effect this has on 
obtaining collision trajectories. 

An important difference in using the planar three-body problem instead of 
the restricted problem is observed when the breakout velocity is determined 
for P3 at L4. As the velocity magnitude ^(0) for t = is gradually increased 
at L4, and P3 moves in the horseshoe regions about the Sun, the location of 
the Earth moves in its orbit, approximately maintaining its 1 A.U. distance 
from the Sun, but shifting its angular position with respect to the Sun. This 
is because as P3 creeps further and further away from L4 approximately 
on Earth's orbit, it can gravitationally perturb the Earth when it moves 
relatively near to the Earth since its mass is now one tenth that of Earth. 
Analogous to the restricted problem, when V{0) gets close to the breakout 
value, the horseshoe region begins to close on itself in a symmetrical way with 
respect to the Earth. However, since the Earth has shifted its location, the 
symmetrical closing point is not on the x-axis as in Figure El for the restricted 
problem, but at another location, the Earth's location, approximately lAU 
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from the Sun. This is illustrated in Figure IT^ where V^(0) = ?/)| = .160 
km/s. a = is assumed. This is fairly close to breakout, for the given 
velocity direction, which occurs for ^(0) = .205 km/s. The final location 
of the Earth when breakout is achieved is in the third quadrant about 20 
degrees from the — y-axis which means that the x-axis is not fixed to the 
Earth but to the mean motion. 
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Figure 19: 7(t), V^(0) = .160 km/s, t G [0, 200]?/ears, x vs y. Sun centered. 

The cases we have examined indicate that the likelihood for collision to 
occur in this problem is similar to that of the restricted problem. However, 
the dynamics of collision is more complicated. 

We now examine a collision trajectory, C12 and the associated dynamics. 
CVl occurs in the breakout state for V^(0) = 205 km/s. This is only discussed 
briefly here as the details can be found in Belbruno & Gott (2004). It starts 
at L4 for t = with a velocity of .205 km/s in the positive x-direction. The 
Earth is located initially on the negative x-axis at lAU distance from the 
Sun, at the origin. 

C12 is plotted in Figure 1201 

In Figure 1201 in the rotating coordinate system P3 escapes L4 moves 
down toward the Earth, then turns around and moves in a retrograde motion 
(clockwise) about the Sun, continuing to the third quadrant, and then in a 
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Figure 20: CL2 from L4 at t = to earth collision at t = 108.628235 years 
in the third quadrant. Sun centered, x vs y. 

posigrade motion (counterclockwise) with respect to the Sun, it circles the 
Sun, moves past the negative x-axis, to Earth collision. Unlike CI = Cll, 
C12 is a posigrade collision orbit. 

Now, as P3 has moved in this collision orbit, the Earth has moved also. 
It has moved in the posigrade direction and then, toward the end in a com- 
plicated motion in the retrograde direction. This is seen in Figure |^ 

The motion of the Earth about 30 years prior to collision with P3 is 
complicated. The final phase of the Earth's motion much enlarged is shown 
in Figure 1221 It consists of many small loops, one for each year, caused by 
the perturbation of P3 as it gets near to collision. In Figure |221 this looping 
motion is shown in large scale about four years prior to collision with P3. 
The Earth is moving in a retrograde fashion in this figure, starting in the 
lower right and ending in collision near the center of the coordinate system. 
The very end of the Earth's trajectory is actually parabolic in appearance, 
as seen in Figure |^ which is too small to be seen in Figure |221 

In the same time frame as in Figure |22l we show P3 on C12 moving to 
collision with the Earth in Figure |2S1 The small black smudge on the lower 
left is the complicated motion of P2 (shown enlarged in Figure |221) prior to 
collision which relative to the scale of C12 is too small to be clearly seen. 



36 




Figure 22: The motion of the earth(note enlarged scale) about 4 years prior 
collision with P3 on C12 at t = 108.628235years. x vs y. 
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Figure 23: C12 4 years prior collision with the earth at t = 108.628235years. 
X vs y. 

The actual Earth collision with P3 is shown in Figure |^ which shows the 
relative motions of the Earth and P3. We have integrated the motion of P3 
through collision to better show the relative motions. In this figure, C12 is the 
larger parabolic type curve, and the Earth moves in the smaller curve. The 
Earth moves clockwise from the left to the right, and P3 moves clockwise 
from the right to the left. The span of the vertical axis is approximately 
20000 km, so that one half of this distance represents the radii of the Earth 
and P3 added together. This implies that actual physical collision between 
the Earth and P3 occurs in this figure when P3 is near the start of its motion 
on the bottom right. We have verified that this is a near parabolic collision 
as in Cll. 

Three- Dimensional Simulation of Collision in an Anisotropic Thin Planetes- 
imal Disk via Random Walk Encounter Dynamics 

We now consider the main model of this paper for the numerical sim- 
ulation of a collision trajectory by a Mars-sized impactor. So, the planar 
three-body problem just considered is now generalized to three-dimensions 
given by (0)). At t = the Earth has the same initial position position as in 
the planar case, and P3 starts at L4. We now more realistically model the 
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Figure 24: Relative motions of the earth and P3 near colhsion. CL2 moves on 
the larger parabolic curve clockwise from right to left, and the earth moves 
on the smaller parabolic curve clockwise from left to right. Time Duration: 
1.9 hours. X vs y. 

breakout of P3 from L4 by the random walk AV accumulation process. 

We assume that A\^'s are imparted to P3 at random times. We will also 
assume that at these random times separate independent A\^'s are imparted 
to the Earth. As is described in the Appendix, we are assuming a thin 
planetesimal disk, and at the random times, two types of velocity kicks are 
applied to both the Earth and P3. One velocity kick is assumed to be in a 
random direction in the plane, and labeled 5V\\ , and the other is perpendicular 
to the plane randomly either up or down and labeled ^V^. The velocity 
kicks applied to the Earth have a subscript of E, and those relative to P3 
have no subscript. In the Appendix we estimate the magnitudes of these 
velocities. The magnitudes for the velocity kicks for P3 are given by in 
the Appendix, and the magnitudes for the Earth are given by (|18|) in the 
Appendix. In this case breakout occurs when k = 14. After these directions 
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are input, the trajectory is propagated for a random time t which varies 
between and 100 years. The process is then repeated. 

Once the breakout state was achieved in the 14th step, and the propa- 
gation terminated after ti4 = 93.1 years, it was found that by extending it 
another 6.9 years to 100 years, no colhsion occurred. We then went back to 
the beginning of the 14th step. A different set of random values were given 
to the velocity kick directions for both the Earth and P3, keeping the time of 
integration to be from to 100 years. Collision again did not occur. We then 
again went back to the beginning of the 14th step, and again picked random 
values, and collision also did not occur. We then made a third additional 
random trial for the velocity kick directions at the beginning of the 14th step 
and found collision did occur when ti^ ~ 4.56 years. 

This yields a probability of collision V = 1/4 (because in one of four 
random trials we succeeded) we discuss this further in the Appendix. (As in 
the discussion following equation IHl our number of trials is sufficient to give 
a rough estimate of this probability which is all we require, and additional 
trials could establish this number to higher accuracy.) 

The collision trajectory 7(t) is plotted in Figure 031 from the initial value 
in the breakout state at the 14th step using the randomly chosen velocity 
kick directions in the final attempt. The time from its initial condition is 
4.5615948 years, and this final portion of the trajectory 7(t) is shown. In 
Figure 1221 it is the upper curve, and the motion is in the downward direction. 
The smaller lower curve shows the motion of the Earth which moves in the 
upward direction. The collision is seen to take place near the x— axis(when 
the center of P3 hits the Earth's surface). Actual physical collision occurs 
slightly earlier when the surface of the impactor hit the surface of the Earth. 
If we continue the trajectory 7(t) through the Earth's surface to Earth pe- 
riapsis, the periapsis distance is only approximately 200 km. The time of 
periapsis is 4.5616056 years. See Belbruno & Gott(2004) for more details. 

The initial conditions for the Earth and P3 at the beginning of breakout, 
4.5615948 years prior to Earth collision are explicitly given in 0. This 
breakout state results from the random walk process previously described 
after 13 velocity kicks at times tj, i = 1,2, 13, where = 4.5615948 years. 
The total time T for the motion of P3 to reach this state is T = ^Ititi = 728.2 
years. The position of the Earth at breakout shows that the Earth has 
migrated in its orbit a considerable distance in a posigrade fashion from its 
initial position on the negative x-axis to the first (upper right) quadrant 
approximately along its orbit. The velocity of the Earth has slightly changed 
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to 29.773 km/s from 29.78 km/s due to perturbations of P3. Even though 
colUsion with the Earth is 4.56 years away, the velocity of P3 is 29.698 km/s 
which differs from that of the Earth by only 75 m/s. 

We describe the collision trajectory of P3 with the Earth. The coordinate 
system is the same that we used in describing the motion of CZ2; that is, 
a rotating coordinate system rotating with the mean motion of the Earth 
about the Sun. It is convenient use Jacobi coordinates q = {qx,(ly,(lz) Q = 
{Qx, Qy,Qz), where q is the relative vector of the Earth with respect to the 
Sun, and Q is the vector from the center of mass of the binary pair Pi, P2 to 
P3 (Belbruno 2004). As with the planar three-body problem we considered, 
we use a rotating coordinate system which initially rotates in the plane of 
the Earth about the Sun, and with the mean motion of the Earth about the 
Sun. 
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Figure 25: Collision between Earth(Lower curve) and Impactor(upper curve) 
in the first quadrant. Earth moves in upward (posigrade) direction, Impactor 
moves in downward (retrograde) direction. Collision occurs slightly above the 
X-axis. Time Duration: 4.5615948 years. Planar Projection, x vs y. 

Along the collision trajectory j(t) of P3 the z-variation between the Earth 
and P3, given by Qz — Qz, oscillates between approximately ±2000 km. This 
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is shown in Figure |2S1 
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Figure 26: Variation of — Qz- Time Duration: 4.5615948 years. 

We show the relative motion of the Earth and P3 near colhsion in Figure 
I2ZI In this figure the time duration is only .53 hours. The orbits of the Earth 
and P3 have been continued beyond collision to get a better understanding 
of the dynamics. The vertical axis spans approximately 10000 km, so that 
actual collision of the surface of the impactor with the surface of the Earth 
would occur near the beginning of the trajectory 7(t) in the upper right 
quadrant. This dynamics is analogous to that of C/2 near collision shown in 
Figure El 

As a final comment, it has been verified that including full solar system 
modeling, as described at the end of Section the process of obtaining the 
collision trajectory 7(t) is perturbed by a negligible amount, and a nearby 
collision trajectory can be constructed. This is due to the fact that at break- 
out, the motion of P3 stays close to lAU radial distance from the Sun. It 
is also noted that in our analysis of the motion of P3 after breakout, it per- 
forms several close flybys of the Earth for several hundred years. During this 
time collision is fairly likely to occur. Our analysis has shown that it is most 
likely to occur soon after breakout, which for the trajectory 7(t) was only 4.6 
years. It has been found that after breakout has occurred P3 generally makes 
very close flybys of the Earth for up to approximately 500 years, and then 
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Figure 27: Relative motions of the Earth (smaller parabolic curve) and Im- 
pactor(larger parabolic curve) near collision. Earth moves counterclockwise 
from left to right, Impactor moves counterclockwise from right to left. Time 
Duration: .53 hours. Planar Projection, x vs y. 

the flyby distance becomes steadily larger. The repetitive flybys increase the 
semi-major axis and eccentricity of the trajectory making collision less likely 
as P3 is no longer restricted to staying as close to lAU radial distance from 
the Sun as it did immediately after breakout. Also, P3 will acquire larger 
deviations in the out of plane directions. Soon after breakout, these effects 
are significantly less pronounced. In this sense, soon after breakout P3 has 
a good chance to creep into an Earth collision; however, if it does not have 
one relatively soon then the flybys themselves will eventually make collision 
less hkely. Overall, as we have shown, the probabihty of coUision with Earth 
rather promptly after breakout is appreciable (of order 1/4). As discussed 
previously, running additional cases could establish this number to greater 
precision. 
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4 Discussion 



We have shown that the giant impactor could have formed at L^iL^) and 
then escaped on a creeping chaotic trajectory to impact the Earth, with a 
near paraboUc encounter in agreement with simulations. 

We note that there are difficulties if the giant impactor came from a 
location other than L4(or L5). To illustrate this assume that it came from 
elsewhere. 

Since the Earth's orbit and Venus' are nearly circular and co-planar even 
at the current epoch, after 4.5 billion years of perturbations, this suggests 
that the early disk of planetesimals in the neighborhood of the Earth was 
quite thin and that the planetesimals in the disk were in orbits that had low 
eccentricity (e) and inclination (i), e ~ i <^ 1. The critical impact parameter 
for coUision with the Earth for a small planetesimal is 

bm = rEii + [vi/v;jy/\ 

where Vpec is the peculiar velocity of the planetesimal, i.e., Vpec ~ e x Vorb ~ 
e X 30 km s~^, and Ves is the escape velocity from the surface of the Earth 
(see Appendix). If V^ec/Krb < 0.004, then bm ~ ?^E(Ks/V^ec) > e x lAU ~ 
i X lAU ~ {Vpec/Vorb)'i-AU, and the planetesimals whose semi-major axes 
are within a distance of bm of the Earth's distance from the Sun of lAU will 
likely suffer collision with the Earth within a short number of years since we 
expect the orbits to be chaotic, and the impact parameter with the Earth is 
less than the critical impact parameter bm for collision with the Earth. This 
will clear out a region of ±6^ around lAU Except for planetesimals in stable 
orbits around L^forLr,). Planetesimals at nearly lAU from the Sun-and not 
at L4(orL5)-will quickly be accreted by the Earth, before having the chance 
to grow large by accretion themselves. T.R. Cowley (lecture at Univ. of 
Michigan 2002) has noted this problem, saying that "Advocates of the Big 
Whack hypothesis usually say that the impactor must have been formed near 
the Earth. This is neither probable nor impossible. It is not probable because 
the Earth could have readily swept up materials that would have formed 
the other body. It is not impossible because we do not know the precise 
conditions of the accumulation of the Earth, and cannot say how improbable 
assembly of the putative impactor near one astronomical unit really was.") 
The stable location at L^^oi L5) answers the probability question, offering a 
reasonably likely scenario for forming the giant impactor near lAU without 
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the material first being swept up by the Earth. Once this cleared out region 
of ±6m has been established, there will be no further quick accretion onto 
the Earth, because the planetesimal's orbits will not take them to within an 
impact distance bm from the Earth. Then they will have to diffuse in by 
two-body relaxation-from perturbations by other planetesimals and planets. 
This two-body relaxation process will slowly put planetesimals into the gap 
region again and there will be be quick accretion from the gap. The giant 
impactor is expected to be one of the later impactors to hit the Earth because 
the successful simulations of the formation of the Moon start with the Earth 
already at nearly its current mass, showing that its subsequent accretion 
(after the giant impactor hit) is assumed to be small (Canup & Asphaug 
2001). The giant impactor should also be expected to be one of the latter 
impacts because planetesimals, including the Earth, grow by accretion with 
time and that would have also allowed more time for the giant impactor to 
have grown by accretion itself. 

If the giant impactor is one of the latter impactors as argued by Canup & 
Asphaug (2001) (after most accretion for the Earth has been completed) then 
if it is not from L4(or L5) it must originally come from either significantly 
outside lAU or significantly inside lAU. But then it would violate one of 
the key advantages of the great impactor theory: namely, point (3) in the 
Introduction, which explains why the Earth and Moon have the same oxygen 
isotope abundance- namely that the Earth and the giant impactor came from 
the same radius in the solar nebula. Meteorites from different neighborhoods 
in the solar nebula(those associated with parent bodies of Mars and Vesta 
for example) have different oxygen isotope abundances. The impactor theory 
is able to explain the otherwise paradoxical similarity between the oxygen 
isotope abundance in the Earth and the Moon combined with the difference 
in iron. 

The Earth has oxygen isotope abundances that are an average over all the 
planetesimals it has accreted-some initially from inside lAU and some from 
outside. A giant impactor forming outside lAU and drawn in by two-body 
interactions would have oxygen isotope abundances intermediate between 
Earth and Mars and therefore not identical with the Earth. Standard giant 
impact theory has the Moon formed primarily out of mantle material from 
the giant impactor material (see (Canup 2004)). The rest of the material in 
the giant impactor is absorbed by the Earth and the iron core of the giant 
impactor eventually finds its way into the Earth's core, leaving the Moon iron 
depleted relative to the Earth. The Moon has been found to have a small 
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core and this is assumed to be from giant impactor material. (See (Wanke 
1999) for a discussion of how the giant impactor theory can accommodate 
this.) If the Moon derives from giant impactor material then it would have, 
the theory proposes, isotopic abundances identical with the Earth if it was 
formed near lAU and this is observed to be the case (Clayton & Mayeda 1996; 
Wiechert et al 2001; Wanke 1999; Lodders & Fegley 1997 . But if the giant 
impactor came from significantly outside outside lAU its isotopic abundances 
would be significantly different from that of proto-Earth. Furthermore, since 
Ml is only 10% the mass of the Earth, this would pollute the proto-Earth's 
isotopic abundances with only a 10% contribution from the giant impactor. 
This would give the Earth and the Moon different isotopic abundances, if the 
giant impactor came from significantly outside lAU. A similar trouble occurs 
if the giant impactor originated significantly inside lAU, if the oxygen isotope 
abundances inside lAU are heterogeneous as well. (At present we have no 
meteorites in our possession whose parent bodies are thought to be Mercury 
or Venus. So we currently have no data for oxygen abundances inside lAU.) 

On the other hand, consider what happens if the giant impactor originated 
at L4(or L5). It is in a stable orbit, so it is not immediately accreted onto the 
Earth, and can grow large and hit the Earth later, alleviating the problem 
mentioned by Cowley. It sits nicely at lAU and accretes exactly the same 
type of material the Earth does, some diffusing from outside lAU, and some 
from inside. The integral of the oxygen isotope abundances of the accretion 
should be identical with that of the Earth. Eventually, perturbations kick 
the giant impactor out of its stable orbit and it collides quickly with the 
Earth. When the giant impactor hits the Earth and kicks out the Moon, 
since the Earth and giant impactor have identical isotope ratios, the Earth 
and Moon should have identical isotope abundances even though the Earth 
and Moon are polluted to different extents by giant impactor material. This 
is an advantage to the giant impactor model, producing automatic agreement 
with proposition (3) of the giant impactor model. Since this is one of the 
latter accretion events for the Earth in terms of the accumulation of its mass, 
the oxygen isotope abundances for the Earth and Moon will not be further 
significantly changed by post-giant impactor accretion. 

Thus we propose the following scenario. 

Debris remains at L4 (as the Trojan asteroids prove). From this debris a 
giant impactor starts to grow like the Earth through accretion as described 
above. As the forming giant impactor reaches a sufficient mass (~ -ImEarth), 



46 



it gradually moves away from L4 through gravitational encounters with other 
remaining planetcsimals and it randomly walks in peculiar velocity. It grad- 
ually moves farther and farther from L4 approximately on the Earth's orbit 
in a horseshoe orbit at lAU, until it acquires a peculiar velocity of approxi- 
mately 180 m/s. The giant impactor then performs breakout motion where 
it performs a number of cycles about the Sun, repeatably passing near to the 
Earth. In a time span roughly on the order of 100 years it collides with the 
Earth on a near parabolic orbit. 

We present here a mechanism for the origin of a Mars-sized Earth im- 
pactor and describe the path it would take to arrive at Earth collision via 
a special class of slowly moving chaotic collision trajectories. The analysis 
shows that Earth collision along these trajectories is likely. Approaches for 
further work are discussed in the Appendix. 

Note added in proof 

As we have discussed, the giant impactor could have grown up in a stable 
orbit at Earth's L4 (or L5) point where a stable orbit is possible and an ob- 
ject could remain and be able to grow by accretion without hitting the Earth 
early-on. We expect this phenomenon could occur when there was a thin 
disk of planetcsimals (in nearly circular orbits). We have noted that Sat- 
urn's rings are an example of such a thin disk of planetcsimals (in this case, 
chunks of ice plus some dirt) observable today. Saturn's regular icy moons 
(inside the orbit of Titan) are all in nearly circular orbits of low eccentricity 
suggesting that they formed out of a thin disk of planetsimals(ice chunks) 
rather like Saturn's rings today only larger in extent. In such a situation we 
might expect our scenario to operate. Therefore it is quite interesting that we 
can find examples of objects at L4 (or L5), or escaping from L4 (or L5) in the 
Saturn system. Saturn's moon Helene co-orbits at the L5 point (60° ahead) 
of the larger moon Dione. Helene has a largest diameter of 36 km and Dione 
has a diameter of 1120 km. Saturn's moons Telesto (diameter 34 km) and 
Calypso (diameter 34 km) occupy both the L4 and L5 points relative to Sat- 
urn's moon Tethys (diameter 1060 km). We would say that Helene, Telesto 
and Calypso originated in a planetesimal disk(of ice chunks) at these stable 
Lagrange points and have grown in place there surviving till the present with- 
out colliding with Dione or Tethys. The rest of the planetcsimals (ice chunks) 
have accreted onto the regular moons of Saturn. (Saturn's rings themselves 
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lie inside the Roche limit where the formation of large objects is forbidden 
by accretion.) While these Lagrange moons are small relative to the primary, 
growth of larger objects with respect to the primary is also possible. Saturn's 
moons Epimetheus (119 km diameter) and Janus (179 km diameter) co-orbit 
in horseshoe orbits just like the one we found for the giant impactor near 
breakout (Figure 5). We would say that Epimetheus formed at a Lagrange 
point of Janus and grew along with it by accretion from the planetesimal 
disk. Later perturbations by other planetesimals kicked it out into a horse- 
shoe orbit just short of breakout. Thus, an object (Epimetheus) nearly as 
large as the primary (in this case Janus) can form and end up in a horseshoe 
orbit. Just a little more perturbation and Epimetheus would achieve break- 
out and likely collide with Janus. These provide examples of the phenomena 
described in this paper that can be observed today. 

A similar pair of co-orbiting objects in horseshoe orbits in another so- 
lar system could be easily detected using stellar radial velocity data. This 
would appear to be a planet in circular orbit about the star whose mass was 
observed to mysteriously vary. The mass variation would be approximately 
sinusoidally in time with a period significantly longer than the orbital pe- 
riod of the primary. For example, if the secondary had a mass 0.1 times 
that of the primary (like the giant impactor) then this would show up as a 
nearly sinusoidal variation of 10% in the deduced mass of the primary. If 
the two were nearly equal in mass there would be a 100% variation in the 
mass. (When they were near each other at one end of the horseshoe orbit 
the effective mass perturbing the star would be nearly doubled, and when 
they circulated to be on opposite sides of the star their perturbation would 
temporarily vanish.) We should have a look among the known cases of extra- 
solar planets for such cases. Granted, we are currently able to see only gas 
giant planets(which may have even migrated inward) rather the terrestrial 
ones we are considering, but still it would be interesting to look. If one found 
such a case, it would be easy to prove. 

Also of particular interest is the Earth co-orbiting asteroid 2002 74^429 
which is in a horseshoe orbit relative to the Earth. Of course, in addition to 
the giant impactor there can be other Lagrange-point debris particularly at 
the other stable Lagrange point not occupied by the giant impactor. This 
material may have been kicked out early-on by other planetesimals or by 
the giant impactor itself as it escaped into a horseshoe orbit. Does any of 
this material survive to the present day? 2002 AA29 (diameter < 0.1 km) 
is in a horseshoe orbit at lAU, virtually identical to the horseshoe orbits 
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found by us in Figure 5. This asteroid approaches the Earth closely (3.6 
million miles away) once every 95 years while circling the sun at lAU. It was 
near one of these close approaches that it was discovered in 2002. After a 
number of cycles, it is briefly captured for a period of 50 years as a quasi- 
satellite of the Earth, before returning to the 95-year horseshoe orbit cycle. 
Its rather large inclination (10.7°) saves it from collision with the Earth. 
This object may have originated near L4(or L5), and have been kicked out 
into a horseshoe orbit (perhaps by the giant impactor itself). If that is so, 
it could be composed of the same material that also formed the seeds for 
the Earth and the giant impactor. A sample return from this asteroid thus 
offers the possibility of obtaining some primordial material from the same 
reservoir that produced the Earth and the Moon. In this case, it should 
have oxygen isotope abundances similar to those found for the Earth and 
the Moon and an iron abundance similar to that of the Earth. The final 
oxygen isotope abundances and iron abundances of the Earth reflect not 
only their seed material (originally from lAU) but also the integral of the 
abundances accreted later, from material originally inside and outside lAU. 
Thus, any slight difi^erences in oxygen isotope abundances would be helpful 
in illuminating the accretion process. It would of course be very interesting 
to measure the age of a 2002 sample. On the other hand, if the sample 
has oxygen isotope abundances identical with the Earth and the Moon, but 
is poor in iron like the Moon, that would suggest it was part of the splash 
material kicked out by the giant impactor at near escape velocity which did 
not coalesce onto the Moon but rather ended circling the Sun at lAU and then 
became trapped at L4(or L5) where it moved for perhaps a considerable time 
before flnally being kicked out. If the sample has completely different oxygen 
isotope abundances from those of the Earth and Moon, that would indicate 
an origin elsewhere in the solar nebula (not at lAU) and we would then have 
to explain how it somehow got perturbed into a low eccentricity horseshoe 
orbit at lAU. (Most Earth-crossing asteroids perturbed into their current 
orbits from the main belt should have much larger eccentricities according 
to Ipatov and Mather 2002.) Bottke, ct al (1996) have previously suggested 
that low eccentricity objects near lAU could have an origin tracing back to 
the Earth-Moon system, and radar results suggest (Ostro, et al 2003) that 
2002 AA2g has a high albedo which supports this hypothesis (according to 
Connors et al 2004). A sample return from asteroid 2002 AA29 is thus of 
particular scientific interest and may provide important clues as to the origin 
of the Earth and the giant impactor that formed the Moon. 
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Appendix: Planetesimal Dynamics in a Thin Disk and Random 
Encounters with the Great Impactor and Earth 



We estimate the magnitudes of the velocity perturbations on the giant 
impactor and the Earth due to encounters with other small planetesimals 
with a back of the envelope calculation. Such velocity perturbations drive 
the impactor of mass Mi into breakout. 

Let /i be the typical mass of the remaining planetesimals (or more pre- 
cisely the rms mass observed in the distribution) (Note: This definition of /i is 
different than the definition of in the discussion of the three-body problem, 
equation P). Consider the disk of planetesimals to be of thickness 



where are the typical peculiar velocities of the planetesimals. and include 
all radii nearer to lAU than either Venus or Mars. The volume of the disk is 
then 2{y^/VoTh)'^{-^'o){lAUY . Let the planetesimals have peculiar velocities 
of order V^, or equivalently orbital eccentricities e(and inclinations i) of order 
y^l|yorb ~ .0006, where Vorb = 29.86 km/s. The thickness of the disk is then 
of order SOte- For distant encounters, with impact parameter b, the deflection 
of the planetesimal by the impactor is 



where be = 2GM//V"^, and m/ is the mass of the giant impactor. 

These are hyperbolic encounters, and the planetesimal exits the encounter 
with the same magnitude of peculiar velocity but changed in vector direc- 
tion. Thus, for nearby encounters, (5V^/V^ can be at most 2. Momentum is 
conserved so the kick in velocity 6V received by the giant impactor is given 



The collisions are independent so the velocity kicks add in quadrature, 
giving a random walk in velocity space with time. Since the disk is thin 




by 



6V ~ 6V^fi/Mi, 



so for a single collision 
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with V^/Vorb ~ .0006, the half-thickness of the disk is .0006AU - 14r^^, 
most of the time h 3> 14r£; and the encounters are mostly in the plane(i.e. 
< 6_L < IAte, while < 6|| < lAU) and 

implying that the major component of the 5V is parallel to the plane. The 
number of planetesimals is n = Mdisk/l-i' so considering the geometry of the 
disk (seen edge on it is a horizontal strip of thickness 2{V^/Vorb)^AU)), the 
number of collisions with impact parameter b (b > IAte) from that strip 
within a time t is: 

SNr^n [2(V;/K,-6)vr(.85)(lA[/)3]-'4(V;/Kr.6)(l^?7)c^6i;,t 

The relevant area is Av{lAU)db, where v = V^/Vorb-, because the disk has 
thickness of 2^(174^7) and there are two vertical strips of width db and height 
2v{lAU) at a distance b from the giant impactor(one to the left and one to 
the right). The range of the impact parameters we should integrate over is 
approximately to lAU. 

The impact parameter for physical impact is calculated from the param- 
eters of the hyperbolic encounter. When the small planetesimal hits Mi it 
will have a velocity V^, where .hVg — .5V^ + GMj/rj. At the maximum 
impact parameter for physical collision the planetesimal will just graze Mj 
at the perigee of its orbit, so that Vs will be tangential at that point and 
its angular momentum per unit mass will be L = Vgrj — V^bm (or equal to 
what it had initially). Thus, b^ = tjVs/V^ and 

brn ~ rj{l + [2GMj/rjV^]y/' ~ rj{l + [bjrj])'/' ~ KiirVbl] + [rj/b,]y/\ 

and since the escape velocity squared V^^ from the surface of the giant im- 
pactor is much greater than V^, then rj/bc — V^^/V^g <^ 1 and bm ~ {bcriY^'^. 
(We are considering here only gravitational encounters, not physical colli- 
sions, so bm will be the minimum impact parameter for our integration. Of 
course, there will be some contribution to (5V||)^ from direct collisions and 
accretion but this is difficult to calculate and we will ignore this contribution 
in this simple treatment. Since {bcViY^'^ is small relative to lAU, this contri- 
bution to {5VY will be negligible in any case. The effects of direct coUisions 
on (5Vj^)^ may be significant, but again difficult to calculate, so for simplicity 
we are ignoring them and only considering gravitational encounters.) 
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Thus, we will integrate from {bcrjY^'^ to lAU, and in the regime we are 
interested in f ~ .0006, (6cr/)^/^ < lAU, and > lAU, so we are in the 
regime where {bcriY^"^ < b < lAU and: 

rlAU 

(SViif ~ iMd^sk/^^)[2v7^i.85)ilAUf]-'Av{lAU)V^t / [V^fxyM]]db. 

(11) 

This yields. 

Breakout is achieved when (5V||)^ ~ (.006)^V^^j, or after a time 

tbreak ~ 7.65 X lO'^-'' {fi/ Maisk){Mj / fifyvs . 

Since v = 0.0006, then 

hreak ~ 35, 400i/rs(^/Mrfi,fe) (M,//i)2. (12) 

At breakout we might expect M^isk ~ 0.3M/ since we want subsequent 
accretion onto the Earth to be inconsequential relative to Mj. Canup and 
Asphaug's successful simulation has the giant impactor hit a nearly formed 
Earth. For example with /i ~ M//300 we expect n ~ 100 other large plan- 
etesimals of mass ~ M//300 and tbreak ~ 35 million years. Long enough 
for iron cores to form in the Earth and the giant impactor as required and 
in agreement with radioactive halfnium-tungsten chronometer results (see 
Canup (2004)). 

The same calculation as done leading to equation ^2 above could be re- 
peated with Me replacing Mj and one would see that 

(13) 

Thus the rms peculiar velocity acquired by the Earth due to velocity kicks 
in time t is 1/lOth as large as that acquired by the giant impactor. 
Returning to the giant impactor we see that for individual collisions 

mr-mribi/q), 

where 

I-V^{IAU)/V,rb rV^{lAU)/V,rb 

< 6l >~ / x^dx/ / rfx ~ {l/3)v\lAUf. 

Jo Jo 
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Similar to the estimation of (5V||)^ in (lllj). it is found that 

nlAU 

{5V^f ~ [M,,,fc/(0.85/i(lAf/)2)]2V(l/3)t;2(lAf/)2 / [V^ ^? / {h' M])]dh 

which simphfies to 

{6V±f ~ [{Md^sk/{0.85|2)]{2/3)V^tv^[V;^/{b,rJy/^M]]. (14) 
Equations (HD), dHj) imply 

{SV^f/m\)' ~ {l/3y[lAU/{hriY/% (15) 
Now {bcriY^"^ ~ rj{Ves/Vfj), and it can be shown that (fTH|l reduces to 
(<5V^±)V('5Ml)' ~ il/3)iVy/Ves)[lAU/rj]. 

If ~ .ooo6v;,b, 

{5V±Y/iSVi\f ~ 1.9 X 10"^ 

So when the giant impactor achieves breakout 5V|| ~ .006^^6, SV± ~ .000026V^rf)5 
so if we simulate the random walk by applying random kicks 

5Vii/Vorb ~ .001, 5V±/Vorb ~ .0000044 (16) 

to the giant impactor at random times, after approximately 36 kicks breakout 
should be achieved. If we followed this to its conclusion that would mean 
an inclination at breakout for the giant impactor of i ~ SV±/Vorb ~ .000026 
radians ~ 5 sec ~ .6rE/{lAU), which would keep it on an easy collision 
course with the Earth. This also guarantees that the collision will be nearly 
in the plane of the ecliptic (i.e. ~ 0.6r£;/0.006yl?7 ~ 4 x 10~^ or within 

0.25° of the ecliptic). This should, by consideration of the total angular 
momentum, produce an orbit for the Moon approximately in the plane of 
the ecliptic as is observed. Some debris is ejected so the alignment should 
just be approximate -which it is. Interaction of the debris disk with the Earth 
as the Moon forms can also naturally lead to a moderate tilt of the Earth 
relative to the Moon's orbit as explained by Canup & Asphaug (2001). 
The calculation for ((5Vj_)^ can be repeated for the Earth, yielding 

{SV±)l/{SVii)l ~ {l/3){mOQf{Vorb/VEes)[lAU/rE]--4.5 x 10-^ 



54 



where K = 2GMe/V^, since r^/foc = V^/V§^^ = (.018 km s-Vll.l9 km 
s-i)2 = 2.6 X 10~^6c = 16AAU. 

Thus, by the time the giant impactor has achieved breakout the values 
for the pecuhar velocity of the Earth will be: 

{SVi\)E/Vorb-- -0006, ((5yx)i?/Kr6 ~ .0000013. (17) 

Therefore, we will apply velocity kicks of 

{SV\\)E/Vorb-^ .0001, {6V^)E/Vorb-^ -0000002. (18) 

at random times to the Earth and after 36 kicks they should achieve by 
random walk the values given by (jl7j) . As can be seen, the movement of the 
Earth is less than that of the giant impactor, so one may say that although 
the motion of the Earth complicates the situation it is still basically the giant 
impactor that is achieving breakout. The Earth is basically a spectator while 
the giant impactor breaks out of its stable equilibrium at L4. 

Let us calculate the minimum impact parameter bm for collision between 
Ml and Me- This is a two-body problem. Impact occurs when the separation 
of the centers of Mj and Me is p = tj+ve- Mj has a peculiar velocity relative 
to Me of order Vbreak- Imagine shooting Mj at the Earth with relative velocity 
Vbreak and a minimum impact parameter b^ so that is just has a grazing 
collision with the Earth. At that point of collision the relative velocity is Vs 
and the angular momentum per unit mass is Vsp = Vbreakbm, which is the 
initial angular momentum per unit mass. Now, conservation of energy gives 

lv^ = lv,lak + [GiMj + ME)/p] 

and this yields 

6„ = p(l + [(V;ynL.)(M, + ME)rE/{MEp}]y^^ (19) 

where Vbreak = -OOGVorb = .18km s~^, Ves = 11.19 km s^^ = 62.1Vbreak is the 
escape velocity from the surface of the Earth. Now, (M/ + Me) /Me = 1.1, 
and p/te = 1.53r£;. Thus, (fT^ yields 

bm = 1.53r£;(l + [3856.4/1.39])^/^ = 80.6rs = .OOSAAU. 

Now, if the vertical peculiar velocity obtained by the giant impactor upon 
breakout is less than .0034 Vorb, as it is in the thin disk case we are consid- 
ering, then it will not be missing the Earth in the vertical direction. It will 
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have a typical eccentricity of .006 and will therefore have a typical impact 
parameter with respect to the Earth in the plane of order .006 AU. So the 
chance of impact in the thin disk case on the first pass by the Earth is of 
order 

V ~ (6„/.006Af/) ~ .57 

or a substantial probability. In fact we observe V ~ .25 as discussed in 
Step 6 of Section 3. In other words we had one success on hitting the Earth 
on the first pass on four attempts. Our earlier tests with the restricted 
three-body problem in the plane had 4 successes out of 13 attempts for 
V ~ .3, as we discussed at the end of Step 2 in Section 3 and given by 
Q. Interestingly, there we picked four directions +x, +?/, — x, — ?/, where 
the values of Hreafc/Kr?) were respectively, .007, .011, .007, .012. The above 
calculation indicates that roughly V oc (l/Hreafc)) so given the values of 
Vbreak wc would cxpcct V ~ .49 for the ±x directions where Hreafc = .007, 
and, in fact, we had 5 runs with two successes at hitting the Earth within 647 
years so we observe V = .4. By comparison, we would predict V ~ .30 for 
the ±y cases where Hreafc = -0115. In fact, in those cases we had 2 successes 
in hitting the Earth in 647 years out of 8 runs giving V = .25. In both the 
±x and ±y cases the results are comparable with our estimates. As discussed 
previously, these are all rough numbers which could be improved by doing 
additional trials. 

(In the case of a thick disk, the perturbations on the giant impactor would 
be of order i ~ .006 = Vhreak/Vorb and therefore it would have a probability 
of missing in the vertical direction of order {bm/ -006) AU, as well as a similar 
probability of missing in the plane so the chance of impacting the Earth on 
the first pass would be of order V ~ [{bm/ -006) AU]"^ ~ .32, which is smaller 
but still appreciable.) 

As we have derived above in equation IT^ 

hreak ~ 35, AOOyear s{Mi / M^isk) {Mj / fi) 

For M^isk ~ 0.3M/ and /i ~ M7/3OO, this gives a breakout time of ~ 35 
million years. 

Since the planetesimals are in orbits with eccentricities e ~ V^/Krb and 
the accretion impact parameter ba > lAU{V^/Vorb), the only planetesimals 
that can hit the great impactor or the Earth are those at a radius from the Sun 
of lAU -bai < r < lAU + bai and lAU-baE <r < lAU + baE, respectively. 
(Recall that the impact parameter for accretion for a planetesimal onto the 
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giant impactor is ~ ^'^^'^^ ~ ^-^ ^ "^"^"^ 

Earth is fe^^^ ~ r^il + V;2/V;2)i/2 _ 522r^ ~ 4.0 x 10^ km > Vf,/VorblAU = 
8.9 X lO"' km.) But once a band of width 2baE is cleared out there will be no 
further prompt physical collisions with the Earth or the giant impactor. 

So the Earth and the giant impactor will quickly clear out that area, but 
further accretion will await scattering of planetesimals into that region on 
a timescale dictated by two-body relaxation among the planetesimals. The 
planetesimal disk closer to lAU than to Venus or Mars has limits 0.86AU 
to 1.26AU. So, to accrete the entire remaining disk, a timescale is required 
similar to that for a planetesimal to random walk up to an eccentricity of 
order 0.2, or equivalently to acquire a peculiar velocity in the disk of order 
~ (0.2)V^,.6. Let us estimate this accretion timescale. 

Let fi be the typical mass of the remaining planetesimals (or more pre- 
cisely the rms mass observed in the distribution). Let the peculiar velocities 
of the planetesimals be ~ 0.0006Vor6. The disk of the planetesimals is to 
be of thickness 2 x 0.0006(l74C/) and includes all radii nearer to lAU than ei- 
ther Venus or Mars. The volume of the disk is then 2(V^/V^,,5)7r(0.85)(l74[/)^. 
For distant encounters, with impact parameter b, the deflection of the plan- 
etesimal upon passing another planetesimal is 

r 2'/^G^^/V^b iib>b, = 2}I^Gy.lVl 

^^^l^^'-X 1 if6<6e, 

where we take into account the fact that the total mass of the two particle 
system is 2/i, the rms relative velocity between the two particles is 2-'^/^V^, 
and is 1/2 of the total change in relative velocity. 
So for a single collision: 



2 



2G'VV^;'&' if ^ > be, 



VI if b < be. 



The collisions are independent so the velocity kicks add in quadrature, 
giving a random walk in velocity space and time. Since the disk is thin, 
be > iV^/Vorb)lAU, {0 < b± < {V^/Vorb)lAU), while < 6|| < lAU so 
the encounters are mostly in the plane and the major component of SV^ is 
parallel to the plane. The number of planetesimals is Mdisk/l^ so considering 
the geometry of the disk [seen edge on it is a horizontal strip of thickness 



57 



2{V^/Vorb)iAU], the number of collisions with impact parameter b (where 
b > {Vn/Vorb)^AU) from that strip within a time t is: 

5N ~ (M,,,fc//i)[2(V^/V;,,)7r(0.85)(lA?7)=^]-M(V;/V;r6)l^f/2i/V^trf6. 

The range of the impact parameters b we should integrate over is approx- 
imately bm to lAU. The rms impact velocity is Considering the 
relative velocities of the two planetesimals and that physical impact occurs 
when their center-to-center separation is 2r^, the minimal impact parameter 
for the physical impact is 

bm ~ 2r,{l + Gi,/r,V^y/' ^ 2r,{l + {bj2'/y,y/\ 

Now, ~ (/x/Mj) i/=^3380 km, and be ~ 2^/^Gij/Vl^ ~ 1.75xl0^km{ii/Mi){0.0006Vorb/V^f. 
So bc/Vf, ~ 5.18 X 10^(/i/M7)2/3(0.0006Krfe/1^At)^ and if the remaining plan- 
etesimals are large but still smaller than the giant impactor (say n/Mj ~ 
1/300) then bjr^ > 1, 6^ ~ 2^/\b,r^y/\ and lAU > 6^ > > {V^/VorblAU » 

(We are considering here only gravitational encounters, so bm will be the 
minimum impact parameter for our integration. Of course there will be some 
contribution to (5V||)^ from direct collisions which is difficult to calculate and 
we will ignore this contribution in this simple treatment. Since bm is small 
relative to be, this contribution to (5V^)^ will be negligible in any case.) 
Thus, we will integrate from bm = 2^/^(6cr^)^/^ to lAU: 

/■6c rlAU 

i5V\^r ~ (M,,,,//.)[27r(0.85)(lAC/)r'(4)2VV,i{ / V;db+ / 2G'V/^>'rf&}- 

Jbm ''be 

This yields, 

{6V^if^iM,,sk/M0.S5)ilAU)]-\A)2'/'iV,/Vor^^^^^ 
which reduces to, 

{6V^if ^ iMd,sk/M0.85)ilAU)]~h6{V^/Vorb){t/yr)[GiJ.], 

yielding, 

iSV\\)VVorb ~ [Mdisk/i0.85)Msun]lWVorb)it/yr). 
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Thus, the accretion timescale is of order 

ta ~ {0.2y[i0.85)Msun/lQMd^sk]{Vorb/V^)yr, 

which yields ta ~ 11.8 milhon yrs [Mj /M^isk]- 

Now, recaU that the breakout timescale is hreak ~ 35,400 years {Mj / Mdisk){Mi / /i) . 
If we want ta > t^reak then 

^Ji/MI ~ 1/300. 

Thus, we expect that the remaining planetesimals would have had time 

to grow large, but would still be smaller than the giant impactor. If Mdisk ~ 
0.3M/, and ^ ~ M//300, there would be of order 100 large (~500 km in 
radius) planetesimals, and the breakout timescale for the giant impactor 
would be of order 35 million years while the timescale for accretion of the 
remaining large planetesimals would be of order 40 million years. This gives 
enough time for the Earth and the giant impactor to form and for their 
iron cores to sink into their centers (in agreement with the estimate of 10- 
30 million years for core formation from the radioactive halfnium-tungsten 
chronometer (see discussion in Canup (2004)). 

So if we want the breakout to occur on a timescale shorter than the re- 
maining accretion timescale, we would want n/Mj ~ 1/300 or greater. Thus, 
we expect by the time the giant impactor is achieving breakout, there would 
be just a fcw(<100) large planetesimals left dominating the mass distribu- 
tion. Reasonable, since the largest of the remaining planetesimals left at the 
time of breakout would have had time to grow large. Still we would expect 
the giant impactor to be the largest of the remaining planetesimals. 

After the giant impactor hits the Earth, and the Moon is formed from 
the splash debris, then the Earth will continue to accrete the remaining 
planetesimals. The accretion cross section for the Earth is of order as ^ 
7r{bcRE) ~ TrrKVEes/V'/i)^ ~ 7r(62rE)^. In Canup and Asphaugh's model the 
Moon is expected to form at 1.2aroche = 3.5rE (much later it drifts out to its 
current location by tidal interaction). The cross section for bringing an object 
inside the Moon's orbit radius is (t^m ~ ■'i'(3.5r£')^(V^g^/3.5V"^) ~ 3.5aE- 
Of those crossing the Moon's orbit, only 1/3.5 will hit the Earth on that 
pass. Of those crossing the Moon's orbit, at the time they cross they will 
have a velocity of V'^ ~ + V|g^/3.5 ~ ^Jes/3.5, and at this velocity 

the Moon's accretion cross section is ~ ^^m(1 + ['^■^'^Mesf^EeslY^^ ~ 
1.167rr^. The fraction that cross the Moon's orbit that impact the Moon 
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will therefore be / ~ 1.167rr^/47r(3.5r£;)^ ~ 1.7 x lO^''. So the number of 
objects hitting the Earth is larger than the number hitting the Moon by a 
factor of (1/3. 5)/1.7x 10~^ ~ 160. If less than 160 large objects are eventually 
accreted then the average number expected to hit the Moon is less than 1. 
So if the number of large objects accreting is <100, there is an appreciable 
chance that all those large objects will hit the Earth and none will hit the 
Moon. This an important advantage for our model since it does not pollute 
the Moon with any additional iron - leaving it iron poor. Indeed, this reason 
is cited by Canup & Asphaugh (2001) in arguing that accretion on the Earth 
and Moon after the great impact be small. The Moon may be expected not 
to gain appreciable additional material. Still smaller planetesimals, which 
make an insignificant contribution to the total remaining disk mass, may fall 
on the Earth and Moon-creating impact sites like Mare Imprium without 
adding significantly to the mass. Since our scenario depends on the fact that 
there will be some debris left in Earth's neighborhood at the time of the 
formation of the Moon by the giant impactor, it is a plus that the Moon 
shows some signs of late impacts itself. 

The parameters here can be considered as a toy model at best. We have 
ignored dynamical friction. Dynamical friction could slow breakout by slow- 
ing the accumulation of peculiar velocity of a massive body, but since L4 is at 
a peak of the effective potential dynamical friction might even speed break- 
out. We have also ignored the effects of momentum transfer perpendicular to 
the plane in planetesimals that actually hit the giant impactor(considering 
only the momentum transfer in the more frequent distant encounters), be- 
cause the latter is difficult to calculate. We have considered pecuhar and 
parallel velocities as if they occurred in a slab ignoring the Keplerian motion 
around the Sun. We need not be married to these particular parameters, 
as they just form a jumping off point. They give us a guess as to the ratio 
of perpendicular velocities to velocities in the plane that might be acquired 
by gravitational perturbations in a random walk scenario. The random walk 
scenario is one that should occur under very general circumstances. Other 
model parameters and assumptions might lead to varying scenarios, but those 
we have shown are a starting point for discussion. A natural continuation 
would be N-body experiments. More complicated mass distributions and 
various eccentricity models for the planetesimals could be considered. There 
is a lot of parameter space to be explored. More elaborate simulations with 
millions of particles including treatment of physical collisions could simulate 
the formation of the Earth and the giant impactor and their growth in a cold 
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disk scenario. After gaseous dissipation was finished and only planetesimals 
were left, we expect some debris to remain at L4 and L5 (like the Trojan 
asteroids) because planetesimals either in, or perturbed, stable orbits about 
about L4 and L5 would stay there. A large object like a giant impactor can 
grow by accretion at L4(or L5). It is a matter of survival: an object in a 
stable orbit at L4(or L5) will survive-not hitting the Earth-and by surviv- 
ing can have more time to accrete other planetesimals and grow large itself. 
One could see how often the second largest object growing near lAU in fact 
started in the Lagrange point debris of the Earth. And finally, in the cases 
where a giant impactor of the type required to form the Moon did indeed hit 
the Earth causing formation of a moon like ours, one could see how often that 
impactor did in fact originate in Lagrange point debris. In other words, what 
is the probability that the giant impactor originated at L4, L5 given that a 
moon like ours(with material of identical oxygen abundance from lAU) is 
formed by collision? 

There are several races going on. The Earth starts forming by accretion 
and as it grows in mass, a stable Lagrange point at L4(and L5) forms. As we 
have shown, as the Earth grows, the region of stable orbits around L4 grows 
in size, so a planetesimal trapped there in a stable orbit would stay there as 
a proto-Earth grew. It would grow larger itself by surviving and accreting 
other planetesimals. It must grow to a mass of order O.IMe, by the time it 
is perturbed out of its stable orbit and achieves breakout. Breakout must 
be achieved before all the remaining planetesimals have accreted onto the 
Earth. After breakout, a collision with Earth on a near parabolic trajectory 
is likely. 

Our paper points out the possibility that the giant impactor that formed 
the Moon could have originated at L4, survived there long enough to grow 
large by accretion, and eventually been perturbed by other planetesimals onto 
a collision course with the Earth. Awaiting numerical simulations capable 
of showing this occurring in detail, we have used examples within our own 
solar system to support our model (see note added). 
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